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ACT 
etric   si-  .     .        ti   ..   is   the   process  of  maintaining  a   set 
.  :fe:  iir-  i  r  of  disturbances.      Geometric 

stabiii:     ti    n    syst  are  :ontrol   of  guided   missiles. 

vc  :.  •'  of  geometric    ~x  ization   use    sin  'n  3   degre  •      f   free 

tegrai        ;  gyros  to   e:  '         Lsh   ir.erti-.l   refer       ■-.  dire       '.on. 

:..  unts   the   gyros   en  a    sta bie    platform  which  isolates   the 
m  m  the   missile   motion.      A   second  method   mc  .:.t  s   t'r.e    gyr 

directly  on  the   missile   anJ    Ln   effect  .-a'-:es   the  missile   the    st    I  L 
.    tfcrm.  this    configuration   angular  moti   n   of  the    vehicle   as 

a    platf*  .-.    I       .   r      severe.      This  .  ay  lead    to   significant  errors 
Ln      aintaining   reference   directj    r.s,    "..hlc\    in   turn   leads   to  error 
in  miss  1  Le    orie        *  .    n. 

Errors  li       Lis   orientation  develop   because   of  the  non- 

iut-itivity  of  three   di::-j:.n'  onai   finite    rotations  and   the    interfering 
effects  of  mi       LI      angular  metier      •      ryro  i  erf  >nnance .     An  ^rror   '  n 
orienteti  r.      in   be   determined   tv    ^onaring   the   angular     ostition  of 
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the  missile  as  se;n  i:.  ir.ertial  axes  v.ith  the  angular  posit.  i.or.  as 
seen  in  tody  fixed  axer. .  The  angular  velocity  of  the  missile  as 
seen  in  body  ax^s  can  be  determined  by  solving  the  missile  ;.erfonna;.ce 
equations.   This  angular  velocity  can  be  expressed  in  ir.ertial  axes 
by  means  of  a  standard  transformation.   The  angular  velocity  can  then 
be  integrated  in  both  coordinate  systems  and  the  resulting  angles 
compared.   Integrating  gyro  gimbal  angles  define  the  missile  orient- 
ation for  guidance.   Com;  n  ;c   of  yro  ;imhal  =sn~le  =nd  missile 
angular  positi  n  as  seen  in  body  axes  l^ter  '     r.o*  ;r  error  in 
orientation. 

.  !       ■  iance  equations  of  the  system  am   "  -li:  v  r  and 
therefore  solutions  were  computed  on  a  digital  conyutir  for  particular 
situations.   The  solutions  showed  that  a  systematic  increase  in 
error  developed  with  time  when  the  missile  was  subjected  to  neri:  :ic 
Lr.t  >rfering  torques.   This  effect  existed  reg  rdless  of  the  control 
syste;r.  dynamics.   When  thenmissile  w^s  subjected  *.o  ran 'I    interf  ring 
torques  the  errors  existed  but  iue  to  the  statistical  nature  of  the 
results  no  further  conclusi  r.s  ..ere  ::  •     Lr.ce  only  z   sin^la  rur.  :  -  s 
made.   Further  ir.vestigatior  of  the  system  using  both  random  and 
periodic  interfering  torque')  is  desirable. 
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Subscript  Symbols 

Syr  bol  Definition 

B  Body  -  refe.-  to  components  in  body  axes 

gim  Gyro  gimbal 

I  Inertial  -  refers  to  component  in  inertia!  axes 

IA  Gyro  unit  input  axis 
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OBJECT 

This   investigation  was   conducted   to  deterrrine   the   effects   of 
angular  base  notion  interference   on  a   geometric    stabilizati   n   system 
using  body   fixed   f~yros.      The   primary  concern  was  the   ability  of  the 
system  to  maintain  a   reference   orientation   in  the   presence   of  errors 
produced   by  the   effects   of   interfering  angular  motion   or:    gyro 
performance   and   the   effects  cf  three   dimensional    finite   rotations 
of  the   base. 
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CHAPTER  I 
INTRODUCTION 

Recent  scientific  and  er  gineering  advances  have  added  satellites 
and  ballistic  missiles  to  the  kinds  of  vehicles  used  by  man.  The 
successful  use  of  the:;e  vehicles  requires  a  means  of  cor. trolling  their 
motion  during  a  cart  o£   their  flight  to  keep  them  on  a  previously 
selected  path.   Ballistic  missiles  and  similar  vehicles  ray  operate 
above  the  earth's  atmosphere.  The  extreme  al'.itu>s  require  rocket 
propulsicn.  The  operating  region  of  these  vehicles  and  the  propulsion 
system  used  gives  rise  to  new  problems  in  implementing  the  control  process 

J  ntrcD  implies  a  means  of  constraining  and  directing  the  vehicle  to 
follow  the  desired  path  from  the  departure  point.   It  is  the  functior  of 
the  control  system  to  direct  the  application  of  the  necessary  constraints. 
The  basic  principles  are  applicable  for  the  entir3  class  of  ballistic 
mi s s i  1°  s  incl  J  ' : . ;  sate  1]  '.  *  j  s . 

R  :  jt  pro;  tiled  vehicles  sxpend  their  entire  stor?  of  chemical 
fuel  duri;  :;  the  first  ;  art  cr   the  flight.   ~his  is  follovel  by  a  relat- 
ively lor  _■   tried  of  :  :  ■■  ■  ■•.-<■<)    fli  -h*  duri:  *  ''  ■  remainder  cf  the  traj- 
ectory. Co   trol  is  applied  during  the  .  v.-efei  part  cr     h^  flight  when 
bl   fuel  is  rapidly  being  converted  into  rrechani  '-     t   .••;  .   The  -\   ' ~    '  '    e 

tve  the  altitude  ?nd  velocity  necessary  to  ;  lace  the  vehicle 
or.  the  ;      '  .-  i  trajectory  -t  the  moment  of  *hrust  teirrr  t tion.  The 
velocity  of  t.  -  vehicH   r  :+  he  determined  with  resvect  to  its  ragnitude 
ani  its  orient  tion.   in  this  resent  it  is  cos  :  !<■>  to  ^rsocia.te  a 


vector  with  the  instantaneous  velocity  of  the  vehicle.  The  task  of  the 
control  system  is  to  identify  the  vehicle  velocity  vector  and  maintain 
its  magnitude  and  orientation  in  a  prescribed  manner.  Most  important, 
the  control  system  must  cause  the  vehicle  to  attain  the  predetermined 
value  of  the  velocity  vector  at  thrust  cut-off.  This  means  that  the 
purpose  of  the  control  system  is  to  control  the  magnitude  and  direction 
of  the  velocity  of  the  vehicle  at  thrust  cut-off. 

In  order  to  control  the  orientation  of  the  vehicle's  'wiocity 
vector,  the  control  system  must  (1)  establish  and  maintain  a  reference 
orientation  in  some  usable  form  (2)  provide  a  means  for  comparing  the 
actual  velocity  vector  orientation  with  the  reference  orientation  and 
(3)  provide  a  means  for  altering  the  velocity  vector  to  make  the  dev- 
iation from  the  reference  become  zero.  The  reference  orientation  need 
not  be  constant  during  the  flight.  It  may  coincide  with  the  desired 
flight  path  or  with  any  other  prescribed  path.  It  is  possible  to  dis- 
tinguish between  an  inertial  and  an  instantaneous  reference  orientation. 
The  control  system's  function  is  to  bring  the  velocity  vector  into 
coincidence  with  the  instantaneous  reference  orientation. 

It  is  desireable  from  the  military  standpoint  to  eliminate  any 
dependence  of  the  vehicle  on  external  sources  of  information.  There 
are  systems  that  are  not  self-contained,  requiring  some  of  the  components 
of  the  control  system  to  be  located  on  the  ground  or  elsewhere  outside 
of  the  vehicle.  There  are  other  systems  that  require  external  radiation 
or  external  reference  points  to  establish  the  reference  orientation. 
Inertial  guidance  systems  have  been  constructed  that  are  completely 
independent  of  external  sources  of  information.  '  These  systems  use 

*  Superscript  numerals  refer  to  similarly  numbered  references  in  the 
Bibliography  at  the  end  of  this  paper. 
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measurements  made  vith  respect  to  inert ial  space  to  establish  the 
reference  orientation  and  to  measure  deviations  from  the  reference. 

Since  the  control  system  depends  on  a  uroc-^s  of  relating  the  v3loci1. 
vector  of  the  vehicle  to  the  reference  orientation  (now  assumed  to  be  in 
inertial  space)  any  uncertainty  in  the  representation  of  the  reference, 
or  in  the  information  that  cor-esponds  to  the  representation,  is  e  qui  vale  n 
to  an  uncertainty  in  the  orientation  of  the  velocity  vector.   The  ability 
to  control  the  orientation  of  the  vehicle's  velocity  is  impaired  to  the 
extent  that  the  reference  orientation  is  uncertain.   To  obtain  high 
accuracy  in  controlling  the  vehicle  trajectory,  it  is  therefore  necessary 
to  maintain  the  refer?nce  orientation  with  high  accuracy. 

The  most  elementary  task  of  the  control  system  is  to  maintain  the 

reference  orientation  v:hen  this  orientation  is  fixed  in  space.   Geometric 

(2) 
stabilization     is  the  ability  of  the  system  to  maintain  the  referenca 

orientation  in  the  presence  of  external  sources  of  interference.   The 
ability  to  track  a  changing  reference  orientation  and  to  command  the  veloc: 
vector  to  acquire  a  prescribed  orientation  with  respect  to  the  instant- 
aneous reference  is  the  guidance  function  of  the  control  system.  However 
a  measure  of  the  capability  of  the  control  system  is  its  success  in 
achieving  geometrical  stabilization. 

3uch  a  stabilization  system  has  been  instrumented  using  three 
single  degree  of  freedom  integrating  gyros  mounted  en  a  gimbaled  platform 
with  their  input  axes  along  three  mutually  perpendicular  directions.   A 

line  schematic  diagram  of  a  three  axis  giir.bal  system  used  for  supporting 

(3) 

such  a  platform  is  shown  in  Fig.  1-1.    The  output  signal  of  each  gyro 

is  proportional  to  the  integral  of  the  angular  velocity  of  the  platform 

about  the  respective  axis.  The  outputs  of  these  gyros  actuate 

gimbal  servo  drive  motors.  These  cause  the  platform  to  rotate  so  as  to 

(4) 
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NOTES:     1.  THE  BASE  MOTION  INTERFERENCE  ISOLATION  GIMBAL  SYSTEM  IS  MADE  UP  OF  THE 
OUTER  GIMBAL  SUPPORTS,  THE  OUTER  GIMBAL,  THE  MIDDLE  GIMBAL,  THE  INNER 
GIMBAL.  THE  ASSOCIATED  DRIVE  MOTORS  AND  THE  ASSOCIATED  RESOLVERS. 

2.  THE  ELECTRICAL  FOWER  SUPPLIES,  ELECTRONIC  UNITS,  COMPUTERS,  CONNECTIONS, 
RACKS  AND  OTHER  COMPONENTS  NECESSARY  TO  COMFLETE  AN  INERTIAL  GUIDANCE 
SYSTEM  ARE  NOT  REPRESENTED  IN  THIS  FIGURE. 

3.  THIS  ILLUSTRATION  IS  BASED  ON  FIG.  4  OF  WRIGLEY,  WOODBURY  AND  HOVORK  A  (2) 
AND  FIG.  6  OF  DRAPER  AND  WOODBURY  (4). 


Pig.     1~J     .ine  schematic  diagram  showing  essential  mechanical  elements  of  inertial  guidance  system 
based  on  rotation  of  the  inertial  reference  package  with  the  indicated  vertical. 


null  the  gyro  signals.   At  null  indication  of  the  gyro  the  platform  is 
essentially  restored  to  its  original  orientation  in  inertial  space. 

Figure  1-2  is  a  pictorial  diagram  of  one  such  instrumentation.   Such  a 

(  2) 

system  may  be   considered   as   three,    single   axis   space   integrators 

operating  in  parallel. 

The  essential  features  of  a   single  axis   space   integrator  are   shown 
in  the   block  diagram  of  Figure   1-3.      This  investigation  was  concerned 
only  with  the   geometric   stabilization  problem.      For  the    remainder  of 
this  paper  it  will  be  assumed   that   the   command   signal  is   zero,   and  the 
actuating  input   to  the    system  is   an  interfering  torque  acting  on  the 
controlled  memcer.      The   perforrr?-nce   of  a   system  ;  ossessing  geometric 
stabilization  may  be   described   by  the    following  steps:    (l)    The   inter- 
fering torque  acting  on  the   control] ed  member  causes   an  angular  rotation 
with  respect  to  inertial  space   about  the    gyro  input  axis.      (2)   The 
angular  rotation  causes  a   precession  of  the   gyro  with   respect   to  the 
controlled  member,   producing  a   signal   from  the   signal   generator,   which 
is  proportional  to  the  angle   of  precession.      (3)   The   gyro  signal  is   fed   to 
the   servo  motor  which   them  applies  a   torque   to  rotate    the   controlled 
member  back  to  its   initial  orientation,  nulling  the    gyro  output   si  r.al. 
System  damping  is  provided  by  compensation  networks. 

The  basic   single   axis   system  then,    consists   of  a    gyro   package 
mounted  on  a   controlled  member  and   a   torqu>    generator.      Th-3    r'ference 
coordinate   system  is   fixed  in   the   controlled  member.      Therefore  maintainin, 
the   controlled  member  orientation  rr.aint-iins   the   reference    coordinate    frame 

With  the   reference   coordinate    systen.  established   within  the  missile, 
the  means  of  controlling  the  missile    to   its  desired   orientation  vrithin 
that  frame  may  be   considered.      It  is  through  control  over  the   orient- 
ation  of  the  missile,    that   control   is  ex-;rcioed   over   the   vehicle's 
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velocity  vector  direction.   In  order  to  function,  such  a  control  system 
must  include  a  mean3  of  generating  the  desired  orientation  of  the  vehicle 
in  the  reference  frame,  a  means  of  measuring  --dual  missile  orientation, 
and  a  means  of  changing  it  to  minimize  or  eliminate  any  deviations 
between  the  d2sir:d  and  actual  orientation.   3ince  this  investigation  is 
concerned  with  stabilization,  it  will  :e  assured  that  the  missile  has 
the  ability  of  attaining  the  desired  orientation  and  the  function  of 
interest  in  the  cci.trol  system,  is  that  cf  maintaining  tnat  orientation 
in  the  ;rese:ic--j  ci"  interfering  torques  applied  t<  the  nissile. 

One  scheme  for  producing  control  torques  for  missile  orientation, 
uses  a  thrust  engine  to  produce  restoring  torques  about  the  pitch  ?ni 
yaw  -<x3s  of  the  mis  >Ile  ani  exhaust  i  .: :  les  or1  '  h  ?  :ircumference  of  the 
airframe  *o  ccntrcl  rotations  about  the  roU  -xis.   This  requires  that 
the  thrust  en/in  be  mounted  in  gimbals,  giving  :  t,  two  decrees  of 
rotaticr.al  f r?e  i  :  ;  bout  trie  ^ero  mo-tent  thrust  j ifne  •  nd  r.  system  for 
changing  the  engine  orientation  in  resj  :>nse  to  command  signals.   Figure  1- 
is  a  diagram  cf  such  a  configuration. 

Looking  at  the  combination  cf  the  missile  ccntrcl  system  and  the 
inertial  space  reference  systen ,  ther?  is  a  means  cf  sensing  deviations 
from  djsired  orientation,  a  torque  generator,  and  s  cor  ■  troilei  member  in 
eac!.  3yst3rr.   In  the  inertial  space  reference  system  the  gyro  package 
senses  deviations  of  the  stable  platform  from  its  desired  orientation. 
In  the  missile  control  system,  the  deviation  sensed,  is  that  of  the 
actual  missile  orientation  from  thst  desired,  ;oth  of  which  are  specified 
relative  to  th-3  stable  platform..   The  stable  platform  is  the  controlled 
member  of  the  ir.ertial  space  reference  system,  an!  the  missile  is  the 
controlled  member  of  the  missile  control  system.   The  inertial  space 
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refei*ence   system  uses  servo  motor'?  as  torque   generators   --nd  the  missile 
control  system  uses  gimbaled   thrust  engines   for  that  purpose,      Slimixation 
of  the   functional  duplication  cf  the   combined    systems   can   result   in 
saving  weight,   and   possible   simplification   of  the   over-all    system. 

Since   inertial   space    reference   is  desired   ar.1    it   is   the   missile 
that   is  to  be    controlled,   the    units  most   reasonably  eliminated   are   the 
stable   platform  with  its   servo   drives,    and    the  missile   control   system 
deviation    sensing   devices.      This   leaves  the   gyro   package   as   the  means 
of  sensing  deviations,    the   missile   as   the    controlled  member,    and   the 
missile   control   system  as   the    torque   gernerator.      The   principle    of 
operation   remains  that  of  a    space   integrator,    but   the   operating  environ- 
ment of  the   gyro  package   has   been   changed. 

In   the   stable   pl'tform  system  for  providing   inertial   space   reference 
the   gyres  are    subjected   to   only  the   motions  of   the    rlatform.      "The 
assembly,    including  the    gyros   and    gimbals,    is    subject   to  three    inputs: 
(1)    concnand   signals  (2)    base  motion   (3)    interference   torques.      Geometric 
stabilization  exi  ts  when   the    gyro   sub-assembly  inside   the     ;imbals 
maintains   its   orientation  with    respect  to   inertial    space    inspit-e    cf  at  y 
base   motion   or   ~he    influence   of    lirectly  actinr  ir.t  ;rf  ^  rence   toroues. 
In  a    sense   the   asserrbly  may  be    considered    ^s  a    base   motion   isolation 
system.      In  particular  it  may  be    regarded   as   rn  active    system   to  isolate 
the    gyro    sub-assembly  against   -ir.  njlar  vibration-,   cf   the    support.      In  a 
rocket  propelled  vehicle,   •.-.•here    the   vehicle    r^nerates   intense    vibration 
during  the    propulsion   interval,    the    gimbsl   supported    base  motion   isola- 
tion  system  greatly  reduces   the   effect   of   the   vehicular  vibration. 
This  enhances   the   ability  of  the    inertial    reference   r.ackage   to  maintain 
the   necessary  reference   orientation. 

(ID 


V.'ith  th^  missile  itself  used  as  the  controlled  member,  the  gyros 
are  rigidly  attached  to  it.  In  this  configuration  the  gyros  are  sub- 
jected  to  all  the  motions   of  the   missile. 

For  op^y^tion   outside   the   earth's   atmosphere   where   aarod;-namic 
forces,      sir  turbulence   and  wini   shepr  are    absent,    the   fissile    structure 
is   still   subject   to   torques   caused     y  lac'-'    cf     /.-•   stry   in   the    exhaust 
flov:  of  the    thrust   3ngir.e ,    -  >;r.3nts  caused    I  y   the   char  ~ir\g  distribution 
of  the   propellent,    and    reactions   en   the    frame   wher    the    thrust    *ngii  - 
nmba]    irive  '   ■•       c\        t:  i    t\  rue?    line    of   ')  5       .    "-   •    relative   to   lhe 

frame.      There'    re,    1         ~        lar  vil  rat:  on    ■:  v'r   r\     n*    of   Lhe    ?yros  vhich 
are    r'    '  il;     a.tta< 
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onv:  r    .     ■•  t  ,    '  r"    ' :    -    ":■  :l-r     otior.3   cf  the   rr.issi:  e ,    -•:-    the   mrir  t- 
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orientation   ^f   the   body  is    "'     1  itely  described  by   the   orientation   of 
its   prii   :ip?.l  ?         .      ""h      torcue    generating   system      rovidas   tha   rr.   ar 
for    •••:'    g   - :    -   !.    ~  rr'  !on   cf  the   missile   ?bout    its    -er  '  sr   of    gravity 

-        ai   ■  s    tc  err   r    -•  [  -r  •  "  ,   -j.->v  1<  ■  ;d    by  the     ■;   •     unit.      The    ]yro 
pack;         :      I * ins  three      '      ni    '     ;  •  -es   cf  fr-     lor-    ir '  egrati  .  •    ~  ro  un4  ' 
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proportional  mechanisms  '-rid  the   affects  of    geometric  cross  coupling 
between  the  missile  equations   of  motion   are   neglected,    such   n    systerr 
becomes  n    combination  of  throe    ind    ;  indent    second   ord^r  sysiams   operated 
it:   parallel.      Cf  course,    tl  is   last   statement   is  r.n   oversimp:*  ficaticn, 
and  while    it   ?ids   ir,  visualization,    the    fret   Is,    that    system  components 
do  possess  dynamic   charsct»risti  :s    arid   there    is   also  coupling  between 
axes    ir    any  real   :  rstem.      Pectus?   cf   tl  is,  certain      r-*    rs    in   svsterr. 
orientation   are    da      3       ■:    .  '  ;  "!i    ar      -1 '  scussed    ir    some    j-tail   Ir.ter. 
Tvr-^e    H'." -->:■■-■  ;rr  ->!   '  '     "'       ';.;•''  I  or.  5 

Consid    r  t!       res;    nse   of  the    svsterr.    iescrited   to   -    small  disturbance. 
A   set  of  inertial   axe     n  •  p    oe    :    fir  ;d    "       :    Lr:oid     ."    with   the   principal 
b<  \'j  axes   of  the   missile   before   the   onset   of  the    disturbance.      The   rotion 
of  the  missile   relative   tc   these    inertia  1   axes    :s  tc   be    ieterrined. 

The  equations   of  .'notion   for  the  fissile   oar.   1  -    c    nv   •  i    ntly  written 
only  in  terms  of  body  axes   components,    since    it   is    onI;»    in  body  sxe:    that 
the   elements   cf  the  inertia    dyadic   for  the   rissile   remain   constant.      If 
such  a   set  of  equations  are  written,    they  make   available    upon  their    solution, 
the   components   in  boiy  axes  of  the    angular  velocity  of  the  missile  with 
respect  to  inertial  space.      The    corresponding  components    ir    the   ir.ertial 
axis   system  will,    in  general,    be   different   and   can  be  determined   by 
means  of  the   standard  matrix  transformation  for  retat^:!   axes  systems. 
That  is: 

(2-1}  '"'(I-B^  =     V-(I-B)B      "      A(I-F)I 

Where  the  subscripts  I  and  B  in  the  vector  equation  (2-1)  refer  to  the 
components  of  the  angular  velocity  vector  as  expressed  ir.  inertial  axes 
and  body  axes  respectively.   The  dyadic  A,    »   is  on  expression  of  the 
rotation  matrix  which  is  derived  in  Appendix  A  of  this  paper.   It  is  made 


(16) 


up  of  elements  involving  the  direction  cosines  relating  the  body 
and  inertial  axes  set.   In  general  these  elements  will  -impend  on 
the  order  in  which  rotations  about  the  thr?e  inertial  axis  occur. 
However,  it  is  shown  in  Appendix  A  that  for  rotations  which  are 
small  in  the  sens?  th-^t  the  product  of  two  angles  is  much  less 
than  one  radian,   he  expression  (2-1)  becomes,  in  matrix  form: 
(2-2) 


,V(I-B)„ 

h 

W(I-B)Y 


•(I-B), 


l . 


~A(I-B)7   A(I-R)V 


-A 


(I-B). 


~A(I-B)V  A(T-B)y   T 


"'•(I-B) 


''•(I-B), 


"(I-B), 


'B 


where  the  subscripts  X— ,  X   Y_,  etc.  identify  the  components  of 
angular  velocity  (V:)  and  the  integrals  of  these  components  (A.) 
in  inertial  and  ^cly   axes  respectively.   Tl  e  2xir?ssion  (2-2)  is 
now  independent  of  the  order  of  rotation. 

The  expansion  of  (2-2)  gives  for  the  comp  ::  nts  of  angular 
velocity  alor.g  the  X-inertia]  axis: 
(2-3) 

W(I-B)    =  '"(I-B).  "  A(I-B)7   ':(I-P)V   +  A(T-B)V   '"'(I-B)7 
71         *3         "I       YB         ^1        B 

Integrating  this   expression    gives   the   angla    through  v;hich  the 

body  has  turned  about  the   X-inertial  axis: 


(2-4) 


j  M(l-B)^dt     -    yW(l-B)YBdt     -  /A(l-E)Zi     T;(I-B)Y^ 

/ 


A(T-B)y       W(I-B)7    it 


(1") 


The  Tirst  term  on  the  right  ir.  (2-4)  is  the  angle  through 
which  the  missile  has  turned  about  the  '.'.   body  axis.  Obviously  if 
the  missile  undergoes  some  arbitrary  rotation  r.nd   th  -  system  of 
Fig.  2-1  is  successful  in  returning  it  to  its  initial  orientation 
in  body  axes  then  the  first  tern:  en  the  right  in  (2-4)  goes  to  zero. 
Since  the  gyro  elements  which  i   •r?lc;  the  error  signal  are  mounted 
on  body  axes,  this' is  precisely  what  the  system  atte-rpts  t  ;  dc. 
This  being  the  case,  in  order  ""or  the  net  "ct-ti  n  ■  c-ut  the  ir.ertial 
axis  to  ~o   to  zero,  it  is  required  that  the  sum  of  the  two  remaining 
t*ms  on  the  right  in  (2-4)  -  tc  :ero.  It  would  a: pear  lhat  this 
'.•.•111  net,  in  general,  occur,  and  some  error  in  c-rier.t-tior.  with  respect 
to  inertial  space  will  leveled  ir.  the  system. 

;ff  ;t  demonstrated  above,  is  essentially  gee::  itric  ir 
nature.   It  is  an  example  of  the  noncemmutativity  of  finite  rotations 
and  the  analysis  can  be  approached  ir.  this  light.  The  nf:  :t  bis  been 
observed  in  the  performance  of  gyro  test  turntables  and  an  analysis 

using  this  somewhat  different  a-proach  has  been  made  by  Goodman  and 

(5) 

Robinson. 

The   equation  (2-4)   rriP.y  be.  interpreted    from  a   slightly  different 
viewpoint  which   h~s  the   virtue   of  indicating  the  nature   cf  the   error 
in  orientation  which  may  develop.      Suppose   that  there   are   two  obser- 
srs;    one   in  the   inerti-  1  reference    frame   and   one   in  the   body  reference 
frame.      Each  observer  concentrates   his   attention  en  his   own   X-axis 
and   is  able   to  record   rotation  of  the   body  with  respect  to   inertial 
space   about  that  axis.      In  addition  the   observer  in   body  axes   can 
cause    rotation  of  the   body  with  respect   to   inertial   space   about   bis 
X-axis.      Nov/  if  the   body  is   subjected   to  an  arbitrary  rotation  the 


(IB) 


observer  ir.   inertial  axes   s^cs   n    rotaMon   nbout  h I r»   X-nxis  which  i? 
the   sum  of  the   terms   o::   the   ri.-ht   in   (2-4)   while   the   observer  in 
body  axes   sees  a    rotation  which  is   3<-ual  to  the    first  term  on  the 
right   in   (2-4).      Let   the   body  now  be   brought   to   rest   at   some   tirre    (t.). 
Suppose   now  that  the   observer  ir.   body  axes   causes   the    body  to   rotate 
about  his   X-axis   only,    so  as  to  make   the  net    rotation  as  he   has   obser- 
ved  it  zero.      During  this   rotation   both  observers    see  exactly  the 
same    rotation.      However,   when  th^   observer  ir.   body  axes   has   finally 
driven  the   boiy   so   that  the   n°t   rotation  as  he   sees   it   is   zero  at 
some   time    (t~),    the   observer  in   ir.ertial  axes   still  observes  a   finite 
net   rotation.      In  fact,    the   observer  ir.    inertia!    axis    sees  a   net 
rotation: 


(2-fO 
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The   advantage   of  the  expression   (2-5)    becomes  apparent  when 
the   concept   of  the   control   system  of  Fig.    2-1  brought   into  the 
discussion.      The   body  axes   set,    the   observer,    ?nd   the  means   provided 
for  rotating  the    body  are   exactly  thi3   control  system.      The   time 
period   (0-t.)   may  be   thought  of  as   the   dynamic    Lime   lag  in  the   syste: 
and  then   the   rotation   observed  by  the   observer  in  body  axis   is  the 
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error  signal  in  the  system.   This  error  si  'rial  itself  represents  en 
error  in  orientation  of  the  tody  with  respect  to  inertial  space,  but 
it  is  tr  nsient  in  nature.   The  system  is  av/are  of  it  and  eventually, 
it  will  be  driven  out.   However,  there  is  in  addition  to  this  transient 
error,  the  deviation  expressed  in  (2-5)  which  the  system  is  not  aware 
of  and  hence,  can  do  nothing  about.   The  advantage  of  the  expression 
(2-5)  is  that  it  expresses  the  error  in  orientation  due  to  the  effects 
of  finite  rotations  independently  of  any  transient  deviations  in 
orientation  due  to  system  dynamics  which  may  exist. 

The  deterrr.ir.at ion  of  the  magnitude  of  the  error  in  orientation 
discussed  above,  and  expressed  in  equation  (2-5)  is,  at  least  in  theory, 
straightforward.   The  equations  (2-2)  taken  together  with  the  per- 
formance equations  of  the  system,  form  a  set  of  equations  which  specify 
the  motion  of  the  missile,  with  respect  to  inertial  space,  in  both 
the  inertial  and  body  coordinate  frames.   It  is  only  necessary  to 
solve  these  equations  for  the  appropriate  rotations  and  werforrr.  the 
subtraction  indicated  by  (2-5),  and  the  error  in  orientation  is 
obtained.   Unfortunately,  these  equations  are  non-linear,  so  that  no 
solution  in  closed  form;  for  the  general  case  is  possible.   However, 
for  a  specific  system,  solution  by  numerical  methods  is  possible. 
Therefore,  if  a  system  which  is  more  or  ]e:s  typical  of  these  likely 
to  be  encountered  in  practice  is  specified,  then  useful  information 
concerning  the  magnitude  of  error  likely  to  arise  in  a  practical 
system  under  the  interfering  effects  of  missile  angular  motion  might 
be  obtained.  The  specification  of  such  a  system  and  the  analysis 
indicated  are  the  purpose  of  the  remainder  of  this  paper. 
Integrating  Gyro  Errors 

Before  proceding  with  the  specification  of  a  system  for  analysis 
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one  other  error  in  orientation  which  can  develop  in  a  systen  such  as 
that  shown  in  Fig.  (2-1)  should  be  discussed.  This  is  essentially  an 
instrument  error  occurring  in  the  single  degree  of  freedom  integrating 
gyro.  It  is  of  interest  because  its  effect  on  the  system  performance 
is,  in  some  respects,  remarkably  like  the  effect  due  to  finite  rotat- 
ions discussed  above.  This  analysis  will  attempt  to  distinguish  between 
the  two  effects  and  effect  some  comparison  of  their  relative  magnitudes. 


For  purposes  of  discussion,  the  geometric  effect  already  described 

s  n< 
(7) 


(5) 

will  be  referred  to  as  the  finite  rotation  effect  '   while  this  new 


error  producing  effect  will  be  referred  to  as  kinematic  drift. 

The  origin  of  kinematic  drift  lies  in  the  manner  in  which  a 

single  degree  of  freedom  gyro  performs  its  function.  The  principles 

(2,U,8) 

of  operation  of  these  devices  have  been  covered  in  detail  elsewhere, 

so  that  a  simole  functional  description  of  their  operation  will  suffice. 

Fig.  (2-2)  shows  a  line  schematic  diagram  of  a  single  degree  of 
freedom  integrating  gyro.  As  indicated  in  the  figure,  the  main  elements 
of  the  device  are  the  spinning  gyro  wheel  mounted  in  a  gimbal,  a  vis- 
cous damper  which  acts  between  the  gimbal  and  the  case,  and  a  signal 
generator  which  generates  a  signal  proportional  to  the  deflection  of 
the  gimbal  with  respect  to  the  case.  In  operation  an  angular  velocity 
of  the  case  about  the  input  axis  (Fig.  2-2)  causes  a  torque  to  be 
applied  to  the  gimbal  bearings  which  acts  at  right  angles  to  the  gyro 
spin  vector.  This  torque  will  cause  precession  of  the  spinning  wheel 
and  its  gimbal  about  the  unit  output  axis •  The  viscous  damper  acting 
between  the  gimbal  and  case  opposes  this  motion.  The  resulting  angular 
velocity  of  the  gimbal  with  respect  to  the  case  is  proportional 
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to  the   angular  velocity  of  the   caue   with   resj>ect  to  inertial   space 
about   the    input  axis,      fence   in  a   given  tin. 3    interval  the   gimbal  wil] 
be  displaced  with  respect  to  the   case   by  an  angle  which   is   proportional 
to  the   displacement  of  the   case  with  respect  to  inertial  space   in  the 
sar:e   interval.      The   signal  generator   rroduces  a    signal  proportional 
to  gimbal  displacement   on:5,   therefore   its   output   is  proportional  to 
the  angular  displacement  of  the   gyro  unit  with   respect   to   inertial 
space   about  the    input  axis. 

The  integrating  -yro  wil]  perform  as  described  above  *s  long  as 
the  only  rotation  of  the  unit  is  about  the  Input  axis.  However,  when 
the  gyro  unit  experiences  rotations  about  the  other  two  ax  s  of  Fig. 
(2-2)  there  are  additional  effects  on  th j  gimba]  motion.  These  add- 
itional effects  which  ^re  discussed  u.  the  following  paragraphs,  are 
the    source   of   kine:     J  ic   drift. 

The   precession   of   the    ^yro. element  under  the    influence   of  an 
input  axis   angular  velocity    carries   the  angular  momentum  vector  cf  . 
the   spinning  wheel  with   it.      Thi  s  produces  a    component   of  angular 
momentum   perpendicular  to   the   spin   reference  axis   of  the    gyro  unit. 
An  angular  velocity  of  the   case   about   the    spin   reference   axis  will 
then  also  cause   precession   of  the  gyroscope  wheel.      The   output    signal 
of  the    gyro  will   therefore   be   in  error  by  an  amount   proportional  to 
the   product  of   the    sine   of  the    gimbal  angle   and   the   angular  velocity 
about  the    spin   reference  axis. 
That   is: 
(2-6) 
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where  in  (2-6)  the  error  is  expressed  as  a  deviation  in  output 
signal.  The  constant  k,  is  a  constant  of  proportionality  and  the 
angular  velocity  of  the  case  with  respect  to  inertia  1  space  is,  as 
indicated  by  the  subscript,  about  the  spin  reference  axis  of  the  gyro, 
unit. 

In  addition  to  the  signal  deviation  produce-]  by  the  angular 
velocity  about  the  spin  reference  axis  a  further  deviation  is  intro- 
duced due  to  angular  acceleration  of  the  ca.:e  with  respect  to  inertial 
space  about  the  output  axis.   Angular  motion  of  the  case  is  transmitted 
to  the  gimbal  and  gyro  assembly  through  the  fluid  of  the  viscous 
damper.   The  viscous  damper  ap  lies  a  torque  to  the  gimbal  which 
resists  motion  of  the  gimbal  with  respect  to  the  case.   If  the  gyro 
wheel  is  assumed  non-spinning  so  that  any  torque  due  to  it  does  not 
act,  the  torque  summation  en  the  gimbal  becomes: 
(2-7) 

gxm   gun     (I-cajOA      d   gim 

where  I  .   is  the  moment  of  inertia  of  the  gimbal  and  rvro  wheel 

gun  * 


about   the   output  axis,    A,T        x~.    is   the   anrular  acceleration   of 
1  '      (I-cajOA 


the 


case  with  rest^ect  to  inertial  scace  about  the  output  axis,  and  c,  is 
the  damping  coefficient  of  the  fluid  in  the  viscous  damyer.   (2-n) 
simply  expresses  the  motion  of  the  gimbal  assembly  with  respect  to 
inertial  space  under  the  action  of  the  single  torque  produced  by  the 
viscous  damper.   It  may  be  rearranged  to  give  the  first  order  equation 
in  gimbal  angular  velocity: 
(2-*) 

I  .  I  , 

gun   ..      .         gim 
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From  equation    (2-£)    it  can  be    seen  th-'t  the  effect    )f  a   uniform 
angular  acceleration  of  the   case  with   respect  to   ir.ertial  space   is 
to  produce,    in   steady  state,   a    constant  angular  velocity  of  the   gimbal 
with  respect  to  the   ca?e.      (The   transient   period    of  this  motion  can   he 
neglected   in  the   present  ca~e   because   the   characteristic   tiir.e   of  the 
gyro,    is   extremely  small  compared   to  the   other  characteristic   times 
of  the    system.)      This   sli;    h- 1  .  ■■  n    gimbal  and    case    obviously  appears 
Li    the    '  :tpu1    signa]     rf  the    signal  genera'    r  •      a    constantly   increasing 
sign.?]    (for  u:  ifcrm  accel    ration).      In   terms    of   signal   ~-vl';'.i"n 
where   the    ratio  of  ginba]    '.   >rtia    to  vise   is    coefficient  ana!   the 
signal   generator   .- •    sitivity  are   I  .a;  '-:    *       for,     a    sin  -1-;    constant: 

(2-9) 

dt 
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The    2    u   tions    (2--')    and    (2-?)   may  r.ov:      •      ombir.    !   to    -ive    the 
total    ievi    :  i        cf  i     ;    gyro   signal  duo    to  motion   about    the    s_:  in 
reference    and  th^   output  axes.      The   resulting  express Ic  )   can  be    given 
the   dimensions   of  sngolar  dispia  :  t   by  dividing  1hr;  igh  by  the 

signal    generator   sensitivity.      This  is   iesiraable    sine-    it  is  nec- 
essary tu  relate   this   deviation  t-j   angular    lis^lacem"    *      f  the  missile. 
If  this   is  don?   the    total   gimbal  angle   due    to  motion   ab'cut  all  three 
axes   becomes: 
(2-10) 
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where  the  second  term  on  the  right  is  the  sum  of   (.1-6)  and  (2-9) 
divide:  by  Lhe  -ignal  generator  sensitivity.   If  the  sensitivity  of 
the  integrating  gyro  unit  from  input  -  xis  angular  velocity  to  gimbal 
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angular  velocity  is  taken  is   unity   'hen   (2-10)   may  be   written: 

(2-11) 

D(Sg), 


tot 
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S      /.    .    ,e) 
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The   similarity  between  this   term  and    (2-5)   will   be   roted.      mh9   gimbal 

angle    in  the   case   where    Input  angular  velocity  and    gimbal  angular 

velocit"  are   kett   ecual   (S.      /. .  •  =1)    represents  the   angle 

1SUUI-BIA»   giir.) 
thrc    j:'  which   tin    systen    is  attempting  tc    turn   the   missile   about   the 

gyro   input  axis.      The    second    term  or.   the    right   is  the   angle    in  tody 
axes   through  which  the   body  wo   Id    be    turned    1:.   the   absence   of   gyre 
errors.      Like   the    expression   (2-5)    the   equation   (2-11)    ;llr.inates 
the  effects   of  any  errors   due    to  time   lags    In   the    system  which  wil] 
eventually  be   driv-jn   out.      The   effects   of   system  dynamics  are   hence 
not   included    In  this  expression.      For  a    given  axis    of  the   missile, 
say  the    X-axis,    the    Input   axis   of  the    gyro   controlling  the    axis  and 
the   body  axis   itself  are    parallel.      Therefore    the   us-    of  the   axis 
system  depicted    in  Fig.  (2-2)    wii3    be    discontinued   and   hereafter  the 
gyro  whose    in.  ut  axis   is   parallel   tc    the    X-axis   of   the  missile   will 
be   referred   to  as   the  X-axis    gyro.      Also   since   the    systen    real3y  acts 
to  control  gyro   gin.bal  angle   and   rot  body  angle    the    signs   in  equation 
(2-11)  will   be    reversed.      This    simply  refers   the   error  to   gimbal 
angle   instead   of  body  angle.      Equation   (2-11)   may  then  be  written  for 
the   X-axis   in  form   similar  to   (2-5)   as  follows: 
(2-12) 
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The  expressions  (2-5)  and  (2-12)  are  now  Identical  in  form. 
They  can  be  combined  to  give  the  total  error  in  body  orientation  as 
it  is  represented  by  the  gimbal  angle.   However,  before  expressing 
this  deviation  a  word  of  caution  is  in  order  concerning  the  quantity 

A  .  .   In  succeeding  discission  the  reference  angle  for  the  system 

gim  °  a  * 

will  be  taken  as  the  gyro  eimbal  an^le,  (A  .  )  .   \s  was  mentioned 

t-j        &  »  \  gim 

above  if  the  gyro  sensitivity  from  ineut  axis  anrular  velocity  to 

gimbal  anrular  vel  -city  is  unity,  then  this  angle  is  the  angle 

through  which  the  system  is  currently  attempting  to  turn  the  riissile 

about  the  particular  axis  being  controlled.   In  the  present  case 

this  sensitivity  is  unity  (1C   integrating  gyro) .     Under  these 

circumstances  the  total  deviation  of  the  body  axes  from  the  inertial 

reference  fram  after  the  effects  of  dynamic  lajs  have  been  removed 

becomes  the  combination  of  (2-5)  and  (2-12) : 
(2-13) 

D(A,  -  A  .  )   =  D(AT  -  O  +  D(V  -  A  .  ) 

x  gxnvx      v  I    B'x    x  E    ginrx 

The  expressions  (2-13)  along  with  the  corresponding  expressions  for 

the  other  two  axes  of  the  system  give  the  total  current  error  in 

orientation  in  the  system.   In  short  they  give  expression  to  the 

difference  in  orientation  between  where  the  system  should  be  going 

and  where  it  is  going. 

In  the  preceding  paragraphs  two  sources  of  error  'which  can  develop 

in  a  system  such  as  the  one  shown  in  Fig.  2-1  have  been  discussed  in 
some  detail.  The  errors  have  been  compared,  and  a  total  system  error 

due  to  both  effects  has  been  expressed.   It  has  been  pointed  out  that 

the  equations  of  motion  describing  such  a  system  are  non-linear  and 

hence  any  attempt  to  determine  the  magnitude  of  these  errors  must 
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resort  to  numerical  methods  of  solution.   This  means  that  a  specific 
system  must  be  chosen.   The  choic?  of  such  a  syster.  and  the  solution 
of  the  equations  of  motion  for  th2  system  cr.  the  dirit;-l  computer 
to  det.^rmin^  the  magnitude  of  the  errors  expressed  ir.  C"*-1"),  (2-12), 
and  (2-13)  will  be  taken  up  in  the  remaining  charters  cf  this  paper. 
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CHAFT2?.  Ill 


TH3   3YS1\2K 


The   analysis  proposed    iii   Chapter  II    requires,    as   a    first   stej  ,    the 
choice   of  a    specific    systei:    such   ■       *  !  2    or  ■■>     -,ho'..T.   in  Fig.    2-1.      Tr 
"r_>!   Tal  this  means  that    -   -  •:  ts,   >"hich   are    'hewn   ?s   blccks    :.r. 

Fig.    2-],    must   be    specified    i      J   lat    their   ;  erf  -  -  •   ■  •    jouat:    •      --;- 

1         ->•■•»>-•->-•-    ti  ~ '!  - 1*    I  r  v      *  *       l       ;     *  *  *  L  st    *****  "-  r*  *  "*  '■    c    of*    >**ri  t  *  *    Jvr*  2*7"*^  *~*s 

7h*'  s   v,"i .",   '  oh.'  -■  !    out    1  r.    til 3    ■  T  ~  '  *  '  '       " ■  *"' •      T)  ^ t ~*  f*c-"o    tr"  ■*?    s  *  s  c if*"*  — 

cation   of  highly    :   -    jic.at-    '      yster    iyr.arric::   v.'ill   a  i  ;    ];ttle    tc    the 
invest}   -*:    n.      Cr    the    other  hand,    '*     'as    !esir?3    tr    r.a":e    1         ystem 
at  least   r  ;..v  tiv        f      '       Lie        hieh  are   st-^  in  *zei   in  this 


*»-    **  l*  *     •  '    -   •  •  ■  •  -* J 


iescril  ^J    in   the    following  paragraphs  wi" 


.f  +  pn  •   '      *         ircl-iH 


::.■;.'    a 


im       •    ** 


..•.-•   -V-      -p  !     -V,*- 


re; 


i ne    .^;>ojrr.   .  . c ..n   r. "t'jr-iJ     .    *ci 

(l)    ~:  c   rrissi^e    or  controlled   Te;r.ber. 

(?)    The    lorque    ..-en  .era  ting   r.ysi    r  ,    '.  •    w'  Lch    i;  rnennt   the 

r  schanism  that   applies   r^storir  •  tcr  ue   to   the 

missile   in   response   to   gyro   sigrals. 
(3)   The    gyro  v-?c)'cr3,   which   includes  both   rate   end   integrating 

gyros. 
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Thea.a  wlH  each   bo    taker,   up  in   turn.      The   necessary   perfonrnr.ee 
equations  ana    the  assumptions   required   to  write   tham  will   be    riven. 
In  addition,   Appendix   ?  of  this  paper  presents  the   detail-id    ieriv?tion 
of  these   performance   equations. 
ThQ  Missile 

The  missile    is  sc-sunei   to  be   s    ri'id   tody.      It<?   orientation  in 
inertial   spac  2c1r.pl jtely  described   by  the   orientation   of   =    set  of 

ortl  ajor.rl  axes   fixed    in  t!.-1    body.      These   5 re   taken  as   pri:   rip?]    a 

with   origin   at    * '         Issile    center    Df    gravity.      Fc   loving   the    reneral 

f  1^N 
aerodynamic    convention  the    X-ax \"    '      * -•.•-.-,*>   *    -.   -*  -.-•"]--   ox  se ,    the 

-  f  tha   mis "i---*1  s   tr:  ' 

•  ocoarion.all;    :      r   forred   tc   as   roll,    ritch,    ~nd   yaw  ■  recce  :tively. 

■p  v*  P  2   ,     •      c  -3 1   *■  r  *    t l  ^        "  ■*        " '   ! ' v    *    '  r.°  rt  \c-  ^    '  ~ x    ?  ^e    *  *  "  *  t    "  n  ^  "  ~r     3 — ^ 

In  an   -~4''ja2  missile    the   burn  in  ■   ~^   fuel    causes  a    "*'  a  *  fi  sart 
shift  of  iV  -    :."        r  of    *rsvitv   an'   char,  r   ..   of   the   -:::r'      of   ^r  ;**tia. 
These    affects   — *e    act   c   r.s~  dared   here       Tu_>   cento**   ~^    '"'v'tv  o*"   * '  ,j 
missile    is   ,js?  -    ?d    fi;:-'    and   +  r  ■  missile   moments    of   ir  >"'-'-    are   assumed 


a  on  r 


;+.sr,f 


strtbi"iiz?ti   ■:"     ^^-f:r    5    of   t 
interest,    translation   of  the   missile   is   rot   eons: !  -  ar  --1 .      Tho-    center 
of   gravity    af  the  missile   is   assumed    fix"-'    '.•     !■  »rtial    sp~ca. 

h'ith  the   above   assumptions  the   equrti  o'  rrotion   of  t!  e  ~is3:l 

a  r*a    e»a silv    otts  *  red    bv    siit^pi  n  o   — ,-~  ■>.  *  ._     -, v-,,*-    »■  \  -,    „•,_..;-<_    -^  _.{.„_   _,.r 

gravity.      In  these    equations   the    arplied   torques    are   th-    "estcrir 

tcre;^    (V.   )    a.\K    li>'   v;-   the    svsten    and   *u«    intsrf  ^rir  ?  *o»-n ;(?*.    .) 

which   represent  vase  mctio:    an~ul  =  r  vibrrticn.      Tueae    acuatior.s 


(^> 


A/otE      Boe//  Axes  (*Sj  ^,^a) 

tfre.  prmcipji  JherT/j/ 
exes. 

X"**tiil  A**S  4n«/  80 «V 
4/e$  C«/"<Ve  o>f  Stiff 
of  Prob/em. 


Missile    CM. 


2, 


MtssjLE  Axfj    Systems 
Fie.    2-1 


ff<&~dl      ?-**•** 
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(Suler's   "-  ntlona)   rroy   t'.en  be  v.-ritten: 
(3-D 

^  X  X 

I  W     -   (T      -  I    )'•'  w     =  M  +M. 

y  y  s         x     z   x         rr.  mt 

y  y 

i  ••;     -  (I     -  T    )'.    '•'     =  M         +  M. 

z   z  y  v'    x  v         m  in* 

Z  7 

V/hers   the   tc  merits   of   inertia  T    ,    I    .    T      ^>r^    -~r!  r.cir  le   -    rents   of   in- 

r<"    *"-    -:SS;,-I        Tl  r  T.jU*.    volcjti   s    V  "  sn1    '■'.'       -i-.    jI-p 

c r>nt      -:.  b?  ->-—-•      of  thi       .^.i3s  r  veloci  ty  of  th*    rr»i  -  si]  3   vit .':. 


tcrcue?   at- cut    the    x'.r    ■    h-J.     s  r."  tv?  ~i-    -■  i  >.      ;;  -.    ]_4    shows    a 

diagram  of  the   i  ne "'•. ani:  .t    by  which  this   is   accomplished, 

Tu_>    restoring  torque    in   r   11    is    supplied   by   auxiliary  exhaust 
PQj**5  jnounted    on   the    ct  rcv>  ■  f  •  •;•■■!  '•<■    of  the   missile    on   o*  r:o.~i''   ■    rid3?. 
The   react  lor    of  exhausting   gas    produce      ?    — r  1?    ?bout    the    r'-ll    s>is 
It   is  e<  '    that   i-rovi:  Lon    '..'   included    for  reversing  ths   thrust 

■j  i  rectj  on  cf  tb  ^  ~  '^    v-'haust      o^*  ^   Ti    **cr   ~  -    *  r~  -  ~*  ^  ~  •*•  ~  *>»,->    ci  ~-A-  ■»--*• 

It   is   assumed    i'     the      resent    syster   th--'     :ontrcl    cf   the   roll 
tcroue    ~  "■*.'"'  *  r  n ;-*   svster;    is    c ronort i  ena 3        ~~    •  i     "  "      *  u  ■      —  L .em 
produces   roll    r33torir.^  moments   in   v^.~   cv.r*   4e    yrc    input    ?'  rr.^ls 
a  c  c "  r J  '  -  -    '      tl   '    f  c "  3 "  •  '■  r.  i   n  quo  t  L  en : 
(3-2) 

M        =3         /  \e,         +   S         /  v^ 

Where   the    sei  "itivitie^   are   the    noil  toreu°    ?ener?tir,  ~   o'    tep  ' 
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sensitivities  in  producing  J"c",1   moment   in   response  to  integrating 
gyro  and   rate   gyro  input  voltares   respectively. 

The  method   for  producing  restoring  torques   in  ;  itch  and  yr.w 
is   somewhat  mere   complicated.     The  missile    rocket  mctor  is   gimhaled 
tc    give   it  two  degrees   cf  rotational    freedom  about  the  Y  and   Z  body 
ax.3.      The   rotation  of  the  nxtor  shout  t\  ?s<3   axns   causes  rot*? tier,   of 
t;.     thrust   lin?    "jo  th?t    it      o  Z    ■  _-ir     ^sses   through   the  missile 
center  of  gravity.      Thi:     rotation   of  t!       thrus*    line    proiuc:S    the 
desired  moment   about    the   axis  tc   le    restore!. 

T-.  the    s'    plest   ,":t'.:  Ve  mechar.isrr    for  r   tatii   ;   t     ■   missile 
motor  about  t:   >.    rimbal    axes   '.-.ould    be   s    :   rid    order.      It   v/c  ild    provide 
'icient   torcue    in   respens?   tc    rvr'    si^r.al   tc    overcame   the    inertia 
reaotioi    of  the   -  -tor  an:     rur.bsl.      In   -  "4:* :  -:..    \   -■-*;.•   is   req  :ired 
to  prevent  oscillation. ^  Such    a      yster:    would   produce   retcr    -imta] 
angle    ':.   response    to   gyre    sir:  al.      If  the   engine   thrust    L~   ass  .red 
constant,    the   squatior    for  restoring  torques    for  small  motor    -Iml'l 
angles  may  he  written.      For  : itch   axis: 
(3-3) 


*m  pcs(e.      ,A      )  iai        Jocs(e        ,A    )F  1  e 

lgu*   rrc  J  *  r^u'  rr.  rsu 


m 

y 

2(DR)        M 
pes  m 

W  2 
n 

pes 

W  2 
n 
pes 

Y.'he  re : 


V,  =  Undamped   natural  frecuency  of  the  missile  r.otor 

n 
pes 

and    its   control   system. 

(DR)  =  Damping  ratio  of  the   missile  rotor  and   its   control 

pes 

system. 
F  =    Engine  thrust. 
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Y* 
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1  «  Distance  from  missile  center  of  gravity  to  motor 
pivot  point. 

pcs(e.   ,A   )  «  Integrating  gyro  signal,  to  notor  gimbal 
angle  sensitivity  of  the  pitch  control  system. 

pcs(e    ,A   )  =  Rate  gyro  signal,  to  motor  gimlel  angle 
r    rgu'  nag 

sensitivity  of  the  pitch  control  system. 
Appendix  B  presents  the  details  of  the  mathematics  leading  to  Equation  3-3 
fcr  both  pitch  and  yaw  control. 
Gyro  Package 

The  gyro  package  contains  three  single  degree  of  freedom  intef.re  :ing 
gyros  and  three  single  degree  of  freedom  rate  gyros.   Fig.  3-2  shows 
how  these  gyros  are  mounted  with  respect  to  body  axes.   Input,  spin 
reference  and  output  axes  are  shown  for  three  gyros   mounted  on  a  base 
which  is  fixed  to  body  axes.   For  the  purposes  of  this  discussion  the 
orientations  shown  represent  those  of  both  the  rate  ^nd  integrating 
gyros. 

The  rate  gyros  are  assumed  to  be  proportional  mechanisms.   Angular 
velocity  of  the  missile  about  the  input  axis  produces  a  proportional 
voltage  from  the  gyro  at  all  times.   The  signal  from  these  gyres  is 
then  of  the  form: 

(3-4) 

e     «  S  f..     \\! 
rgu     rgu(w,e;  x 

Where:   S   #Tr  x   _,.   .    .    .     ,   ..  , . 

rgu(V<,e;  =  ine  input  axis  velocity  to  output  voltage 

sensitivity  of  tha  r?-te  gyro  unit. 
The  integrating  gyros  are  assumed  to  be  Tree  of  drift  effects. 
Otherwise,  the  units  are  assumed  to  be  identical  to  the  unit  described 
in  Derivation  Summary  2  of  reference  (2).  Appendix  B  derives  the 
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performance  equ;  tion  for  the   integrating  gyro  ur.it   in  terns  of  the 

gyro  axes   givin  in   Fig.    2-3«      These   are   then  changed   to  give   the 

eouations  in  terms  of  tody  axes   components   for  the   orientation  of 

the   gyros   shown  in  Fig.   3-3.      The   result    for  the   X-boiy  axis   (roll) 

be  c  ome  s : 

(3-5) 

(CT).        e.        +  e.        +  V.'     e.        =  S      ,.  s        W     +   (CT) .        V 

v      'igu     igu  igu         y     i -a  sg(A    .    ,ej  x  'igu     z 


vViAre*  (r"r) 

i7igu  =  Characteristic   time   of  the   integrating   ryro   (App.    ?) 

e.  =  Gyre   outout   signal  voltage, 

igu 

q 

sg(A    .    ,e)    =   Crvro  signnl    z  :.i--*^r   gimbal    sngle   to  voltare 
gurr  " 

sensitivity. 
Similar  equations  are    giver,   for  the    other  two  body   axes    ir.   A;    endix   F. 

The   equations   (?-l)    through   (3-5),    alcr..-'  with   similar  e.-.pressior.s 
to   cover  all  body  axes,    completely    define    the   performance   cf  the    system. 
These  might  be    conbir.el    in   a    straight  forward   manner  to   give   a    single 
performance   equation   for  each  axis   of  the   ::.issile.      However   in 
programming  for  the   digital   computer  it    is  much   sirrvler  to   treat 
the    component  performance   equations  as    given  above.      Therefore   no 
overall   ;  erformvnee   equation   for  the   axes   are  written. 

Appendix  B   gives  a   collection  of  the    performance   equations   for 
all  components   for  all  three   body  exes.      It   is   convenient  to  make 
certain   changes   of  the   variables   in  these   equations   at  this   point. 

The   equations  of  motion  for  the  missile    (3-1)    about  each  axis 
are   divided   through  by  the  moments   of   inertia   of  th^  missile   about 
their  respective   axes.      This   changes   the   units   of  the   equations    from 
those   of  torque   to  those   of  angular  acceleration.      The   ratios  of  restoring 
torque  and   interfering  torque   to  the  moment  of  inertia  may  then  be 
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defined  a?  new  variables   representing  the   angular  acceleration  of 
the   body  due   to  these   torques  acting  alone. 

The  definition  of  annular  acceleration  due   to   restoring  torque 
can  now  be   carried   into  the   equations    for  the   torque    generating  systerr 
by  dividing  these  equations  through  by  the  appropriate  moments  of 
inertia.     This  puts  both  the   body  and   torque   generating  systerr.  equations 
in  terms  of  angles  and   their  derivatives. 

The   gyro  equations  can  be   expressed   in  terns   of   gyro   gimbal  angle 
by  dividing  through  b;     the    signal   generator   sensitivity.      This   change 
also   carries   through  to   the    torque    generating  system,      'hen  all  of  the 
above   changes   have    been  made    the    overall   sensitivity  cf  the    system 
appears  as  a   coefficient   in  the   equations   of  the   torque    generating 
system.      Appendix   B  performs  the    changes   of  variable    indicated   and 
redefines  the   coefficients  of  the   equations.      The    resulting  set  of 
equations  along  with  the    transformation  equations   (2-2)   define   a   complete 
set   for  the   problem  under  analysis.      These   are: 
(3-6) 


•  • 


w     -  (ir)  v:  W     -  A       +  (S       ,,         V  xA   .        +  S       /,.  "  v  ,. 

x  x  y  z         m  rcs(A    .      A   )    gim  rcs(V:  A    ;  V. 

*  x  gim'   m     °     x  '   m       x 

•  ••  •  • 

w    -  (ir)  v;  v:    =  a      +  a. 

y  y  x  2         m  int 

•  •  •  •• 

W     -   (IR)  W  W     -  A        +  A. 
z  z  x  y         m  int 

"  z  z 

(4)  (3) 

A  +   2(DR)        K         A       +  W  2     A        =  W  2     S        ,,         "    x   A    .        + 

m                    pes  n         m           n         m  n          pcs(,A    .    ,A    ;      gun 

y                 r         pes     y           pes     y  pes  r         gim'  m       °     y 

w  2   s     fv  ;•  x  vi 

n         pcs(Y.,A  ;     y 

pes  r         '  m 

C  *  2<DRW„       i3'  ♦  V  2     A       .  V  2     S       „        -    /.        ♦ 

z  "  yes  m  n         m  n         ycs(A    .    ,A    )   gim 

*  z  yes     z  yes  *  g3.ni     m  z 

-  W  2     S        fw  V  ^  w. 

nycs  ycs(W'Air)      Z 
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(CT),      A,         +  A    .        ♦  V/  A    .        -  W     -  (CT).      W 
vv   'igu  gim  gim  y  gim  x  'igu  z 

(CT),      a\        +  A    .        +WA    .        -  W     +   (CT).      W 
v      'igu  gim  gim  x  gim  y  'igu  z 

(CT),     *A    .        +  A   ,        +  W  A  -  W     +   (CT),      W 

'igu  gim  gim  y  gun  z  'igu  x 

z  z  z 

W       -  W     -  A     W     ♦  A     V/ 

Xj  x         Zj  y         yx   z 

K       «  W     -  A     W     +  A     W 

yx       y       *!  z       z!  x 

v;       «  W     -  A     W     +  A     w 
Zj         z         yx  x         Xj  y 

Where:    S        ,.  7   \  ■   Gyro  gimbal  anele   to   restoring  angular  accelyr- 

rcs(A    .    ,A    ;  ^        <= 

gim'   m 

ation  sensitivity  of  the    roll   control  system, 

S        /,.,A    )  =  Angular  velocity  to   restoring  angular  acceler- 

rcs(V. '   m  & 

ation  sensitivity  of  the   roll  cot.trol   system. 
It  will  be   noted   that  the  equations   (2-2)   have   been  modified   by  dropping 
the   subscripts   indicating  body  axes  components.      This   is  in  agreement 
with  the   system  equations   of  motion.      Angles  and  their   derivatives 
which  are    referred   to   the    inertial  axes    s*t  are   identified   by  the 
subscript   (i).      The    sensitivities  defined    for  the   roll  control   system 
above   have   identical    counter  parts   in  the   pitch  and  yaw  control   systems. 

The   constants    required   for  the    set  of  equations    (3-o)   must  be 
specified.      The   complete    3'it  of  these   constants  for  the   system  are 
given  in  Table   III-l.     The   inertia   ratios  ware   arbitrarily  chosen. 
They  represent  what  might   be    found   in  a  missile   of   small  to  medium 
size   such  as  the   second   stage    of  a  two  or  three   stage   vehicle.      The 
integrating  gyro  characteristic    tines  are   those   for  the  MIT  Instrumenta- 
tion  Laboratory    10      -j^rro.      The   natural  frequency  and  damping  ratio   for 
the   torque   generating  systems  were  arbitrarily  chosen.      The   values 
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Constant 

Value 

Units 

(IR)X 

0.19S1 

(IR)Y 

0.8912 

(IR)Z 

-0.925 

S(rcs)(Agija;A) 
S(pcs)(Agim;A) 
S(yc3)(Agira;A) 

-1.55 
-27.4 
-28.0 

/Sec2 
/c     2 

/See2 

S(rcs)(V:-A) 

-  3.11 

/Sec 

S(pcs)(V:;A) 

-  5.5^ 

/Sec 

i 

S(ycs)(W;A) 

-  5.44 

/Sec 

w 

n 

25 

Rad/Sec 

pes 

w 

25 

Rad/Sec 

yes 

(DR) 

pes 

0.5 

(DR) 

yes 

0.5 

(CT). 

.003 

1 

Sec. 

TABLE  3-1 

TABLE  OF  SYSTEM  CONSTANTS 
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chosen  are   representative   of  what  might  be  expected   from  a  hydraulic 
servo  in  combination  with  the   engine   inertia.      The   system  sensitivities 
which  have  all  been  combined   in  the  torque   generating  systems  remained 
to  be   chosen.      These  were   chosen  on  the   basis  of  a   linearized   stability 
analysis  for  the    system  which  is  discur^ed   telovf. 

In  an  effort   to  obtain   sore   information   on    system  stability,    the 
equations   of  moticn   for  one  axis  were   combined.      All   coupling  terms 
from  the   remaining  axes  were   assumed   to  be    zero.      The   root   locus   of 
the   characteristic  equ-ticn  of  the   resulting  fifth   order  system   is 
shown  in  Fig.    3-3.      The   calcjlations   leading  to  this    locus  are   given 
in  Appendix  B. 

The   locus  of  Fig.    3-3   shows   two  oscillatory  modes.     As   the   open 
loop  sensitivity  is   increased   from  zero,   the   pole   pair  associated  with 
the   torque   generating   system  moves   toward   the   rieht  half-plane.      At 
the    same   time   the   pole   pair  at   the   origin   due   to   the    integrating  gyro 
and   the  missile   performance   equation  move    to   '.he    left.      The    system 
open  loop  sensitivity   (3        /,  ■•    *)   was   chosen   so  as   to   give   both 

pCSv,  A     .      ,  A     ) 

gun'  m 
pole  pairs  the  same  closed  loop  damping. 

The  preceding  paragraphs  have  served  to  develop  the  performance 

equations  of  a  simplified  version  of  a  geometrical  stabilization  system. 

It  has  been  pointed  out  previously  that  the  errors  under  investigation 

occur  regardless  of  system  dynamics.   This  fact  can  be  demonstrated  by 

varying  the  dynamics  of  the  system  just  developed  and  solving  the 

problem  for  the  same  interfering  input.  Since  the  problem  was  solved 

on  the  digital  computer  the  variation  in  dynamics  can  be  accomplished 

by  rather  elementary  changes  in  the  computer  program.  Therefore  another 

system  is  proposed.   The  body  equations  of  the  sat  (3-6)  were  retained 

but  the  torque  gene  rating  systems  and  integrating  gyros  w-ere  assumed 
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to  have  no  dynamics.  A  set  of  equations  similar  to  the  set  (3-6)  is 
given  for  this  system  in  Appendix  B.  The  programming  of  these  two 
systems  for  the  digital  computer  is  taken  up  in  the  next  chapter. 
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CHAPTER  IV 


METHOD  OF  SOLUTION 


Nori    linear  differential  equations  are   not   in   general  amenable 

(12) 

to  solution  in  closed  form.     Finding  the  exact  solution  of  a  number 

of  simultaneous,  non  linear  differential  equations  is  even  less  probable. 
However,  numerical  techniques  have  been  developed  for  obtaining  partic- 
ular solutions  of  differential  equations,  either  linear  or  ncn  linear. 
The  solution  so  obtained  will,  as  the  name  suggestsv  be  in  number  form. 
For  solution,  the  equations  must  be  stated  with,  numerical  coefficients. 
The  coefficients  may  be  variable,  but  the  manner  in  which  they  vary 
must  be  known.   Certain  initial  conditions  of  the  variables  are  also 
needed.   The  computation  time  required  to  solve  for  th«  dependent 
variables  over  a  given  range  of  the  independent  variable,  will  depend 
upon  the  numerical  technique  and  the  accuracy  required.   For  improved 
accuracy  the  increment  of  the  independent  variable  between  calculated 
points  must,  in  general  be  decreased,  thus  increasing  computation  time. 
The  r.umerical  method  used  in  solving  the  s^t  of  ncn  linear 

differential  equations  (3-6)  which  are  listed  in  the  previous  chapter 

(13) 
is  the  Runge-Kutta  Method  with  the  Gill  modification.   '  Use  of  this 

method  requires  a  known  value  of  each  variatle  and  the  first  derivative 

of  each,  at  the  start  of  each  integration  step.   If  the  differential 

equation  is  higher  than  first  order  the  integration  process  must  be 
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used  as  many  times  as  the  order  of  the  equstion.   For  example  if 
2 

""p  ■  F(x,  y,  z,  t,  -rr*   -rr  )  the  value  of  the  second  derivative 
dt  aL  a1, 

of  y  is  calculated  at  a  point,  say  t  ,  where  values  of  x,  y,  their 

first  derivatives,  z,  and  t  are  known.  The  integration  method  will 

give  the  value  of  the  first  derivative  of  y  for  t  =  t  -  At  where 

At  is  the  integration  interval.   To  get  the  value  of  y  at  the  point 

t  +  At  the  first  derivative,  whose  value  was  calculated  for  t  =  t 
o  '  o 

from  the  previous  time  step  which  started  at  t  -  At,  is  integrated 
using  the  same  process.   The  equations  which  indicate  the  need  for 
solving  for  lower  order  derivatives  are  called  auxiliary  equations. 
The  eight  second  order  differential  equations  mathematically  describing 
the  system,  therefore  require  that  the  integration  n.ethod  be  applied 
sixteen  tires  at  each  time  step.   The  transfer  to  inertial  axes  of 
rotation  rates  expressed  in  body  axes  give  three  more  derivatives 
that  must  ve  calculi ted  and  integrated  for  each  time  step.   The  equations, 
in  computation  form,  along  with  their  auxiliary  equations,  are  shown 
in  Appendix  B. 

The  basic  tool  used  in  solving  the  differential  equations,  given 

in  Appendix  B  was  the  I.  B.  K.  650  mag.-.etic  drum  data-processing 

(14) 

macnine.     The  unit  used  is  operated  by  the  Mathematics  Group  in  the 

M.I.T.  Instru-.entation  Laboratory.   The  machine  used  has  an  internal 
storage  of  two  thousand  re  isters,  two  tape  ur.it s,  floating  point 
arithmetic,  and  index  accumulators. 

The  problem  was  programmed  using  a  floating  address  assembly 

(15) 

routine,  called  Flad,    developed  by  the  M.I.T.  Instrumentation  Lab 

Math  Group.  The  basic  advantages  in  using  this  assembly  routine  are 
the  simplifications  in  coding  and  the  optmizaticn  of  the  program  to 
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reduce   computation  time.      The   differential  equation   routine   was 
available   as  a    subroutine    requiring  only  entry  and  exit   programming. 

The   program  could  be   broken   do'..rn   into  functional  parts  as   follows: 

1.  Obtain  the    forcing  function,   A/.       \,    for  each  axis. 

2.  Compute   the   values  of  the   derivatives  and   store    them  in 
the   proper  registers   tc   enter  the    iiffercntial  equation   subroutine, 

3.  Entry  and    exit  instructions   needed  v.ith  the    subroutine. 

4.  iut  the    required   numbers   in   the    desired   form  and    read 
them  o,.t  of  the   c.mputer. 

5.  Start  the    cycle    over  again  at   the   next   time    step. 

In  running  the   :  roclerr.  the    time   step  needed   to   give   the    required 
accuracy   of   solution   prove  :   to   be    ^nforta^iatJy   small.      The    root  locus 
of  the   linearized    system,    given   in  the    preceeding  chapter  indicated 
thet  the   postulated  system   should   be    stable.      In  wringing  cut   the 
program  a   sinusodal  :'or;ir.g  function  of   15  radians   per   sec::;'}    •  ■"-s  used. 
In  a   linear   system  time    increments    of  one   twentieth  of  a   period   usually 
give  accurate    results.      For  such   a   time    step   the   solution   of  the    system 
equations  diverged.      Using  one  millisecond   time    steps   shewed   that  the 
solution  was   stable   and    some    of  the   variables  changed   rapidly   requiring 
that  the   time   ste;    be   small  for  the   digital   solution   tc   follow  the 
changes.      By  trial  and    err^r  it  vas    found   that   for  a    time   step  of  five 
milliseconds   the   ccrr.uted   values   of  the   variables   agreed   with   the    same 
computed   values  when   one  nilliseccr.d  was   used.      At  ten  billisecond 
time   steps  the    solution   diverged.      Therefore    five   millisecond   time 
steps  were   chosen.      This  time    step   resulted   in   approximatly  twenty  one 
hours   of  computing  time   being  required   tc    simulate   one  minute   of  problem 
time,   using   square   pulse   forcing  functions. 
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Suggested  methods  of  forcing  the  equations  were  by  the  use  of 
actual  acceleration  data,  summing  sine  waves,  and  using  pulses.   Use 
of  data  would  require  reading  data  cards  at  each  computation  point. 
This  would  have  required  the  punching  of  at  least  four-  h  mired  cards 
for  each  second  of  problem  time.   Evaluating  a  sine  t-ikes  150  milli- 
seconds.  For  a  single  sin-;  wave  per  axis  almost  a  half  a  second  would 
therefor?  be  required  each  time  the  forcing  function  was  calculated. 
For  a  square  pulse  only  additions  and  subtractions  are  required  and 
little  time  is  spent  in  the  forcing  function  part  of  the  program.   The 
selection  of  square  pulses  used  in  t1"-!?  ir.vee tig?. tier,  '-'as  primarily 
determined  by  the  computation  time  required. 

The  system  was  programmed  to  determine  the  response  to  three 
forcing  functions. 

1.  An  initial  displacement  of  the  r.-.issile  about  all  three 
axes. 

2.  Feriodic  rectangular  pulses  of  angular  acceleration 
applied  to  a?1  three  ax  -s. 

3.  Random  rectangular  pulses  of  angjlar  acceleration  acting 
on  all  three  axes. 

The  initial  displacement  of  the  missile  about  all  three  axes  was 
used  to  determine  the  transient  behavior  of  the  system.   A  similar 
initial  angular  displace  tent  condition  could  be  realized  in  an  actual 
missile,   tor  a  multistage  missile  the  sejaration  process  would  be 
expected  to  cause  some  rotation  of  the  next  stage.   When  the  propulsion 
system  fired  the  control  system  woul  1  see  the  displacement  as  an  initial 
error  in  orientation. 

Rectangular  pulse3  were  used  for  periodic  forcing  of  the  system 
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primarily  due  to  the  speed  with  which  they  could  be  generated,  as 
noted  earlier.   Use  of  these  pulses  added  no  appreciable  commutation 
time  over  that  required  for  the  initial  displacement  case,   The:;e 
pulses  had  a  magnitude  of  four  tenths  radian  per  second  squared. 
For  each  axis  the  puls^  duration  wa.^  a  tenth  of   a  second  with  two 
tenths  of  a  seconl  between  the  -n'  cf  on-  pulse  ;«n1  the  sta.rt  ;>r 
?.n  "> '  ,h - c .   The  pulses  were  phased  so  only  one  axis  at  a  time  was  being 
forced,   '..her  the  pulse  was  taken  of:  of  one  axi^  and  applied  to 
another,  its  sign  was  changed.   This  gave  the  period  of  the  forcing 
function's  fundamental  frequency  as  six  tenths  of  a  second. 

The  random  pulses  of  angular  accelerati' n  ere  limited  in 
magnitude  fro:;;  zero  to  ninety  mill i radians  per  second  squared,  in 
increments  of  ten.   Pulse  duration  was  limited  to  five  hundred  billi- 
seconds in  increments  of  five  milliseconds.   The  sign,  magnitude  and 
duration  of  each  pulse  was  determined  usin~  digits  from  registers 
containing  calculated  values  of  scn.e  variable.   These  registers 
contained  values  in  floating  point  fori:..   No  digits  less  than  two 
digits  to  the  right  of  th?  first  s' gr.lf  icar.t  figure  v;-re  used.   These 
digits  varied  at  every  tim  2  step.   The  digits  used  to  determine  the 
forcing  function  for  one  axis  came  from  variables  primarily  influenced 
by  motions  about  the  other  axas.   A  single  digit  determined  the  pjlse 
magnitude.   To  determine  the  sign  of  the  pulse,  five  v/as  subtracted 
from  a  single  digit  and  the  sign  of  the  difference  was  used.   The  length 
of  the  pulse  ir.  five  millisecond  increments  v%as  determined  using  two 
digits.   These  procedures  give  a  forcing  function  on  each  axis  which 
is  independent  of  the  mode  of  forcing  on  the  other  two  axes.   The 
couplirg  between  axes  given  by  2uler's  equations  and  the  gyro  equations 
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does  not  however  make  the  angular  motion  -bout  any  cne  axis  independent 
of  the  motions  about  the  other  two  axes. 

Programming  for  the  further  simplified  systems  which  were  mentioned 
in  Chapter  III  was  accomplished  by  removing  parts  of  the  original 
program. 
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CHAPTER  V 

Discussior:  of  results 

This   chapter  discusses   the    results  obtair.ed   from  the    solution   of 

the   system  equations  on   the   digital  computer.      The   discussion  is  concerned 

primarily  with   the   angular  differences  developed    in  Chapter  II.      These 

differences   represent  error?    ir.   indication  of  body  orientation.      The 

difference    cetween   inertial  angle   arid   body  angle    is   the  error   in 

orientation  due    to   finite    rotations.      The   difference   b^tv.eer.   body  angle 

and   gyro   gin.bal  angle    is   the   error  in  orients  ticn   caused   by  gyro 

characteristics.      The    3^.  of   these   differences   (A_-A    .    )    is   the   net 

I      gun 

error  in  orientation  due   to  these   effects.      These   difference-,  do  not 
include   the   deviation   in  orientation   resulting   frorr:  dynamic    lags   in 
the    system. 
Transient  Response 

It  was  desireable   that   sor^e   examination   of   systen    stability  be  made. 
In  order  to    do   this  the    response    of  the    system  to  an   initial    angular 
displacement  was  computed.      No   interfering  angular  accelerations  were 
ap.Jied  during  this   response. 

Figs.    5-1,    5-2,   and   5-3   show  the    response    of  the    system   in  restoring 
from  an  initial  angular  displacement.      In  these    figures   the   gimbal  angle 
is  taken  as   the   reference   variable.      This  was  done   because   the    gimbal 
angle   is  the  variable  which  the    system  attempts   to  null.      The    figures 
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for  the   pitch  and  yaw  axes   show  a   damped   oscillatory  respor.se.      The 
response    in   roll  was  only  run  until  the   gimbal  angler,  on  the   other 
two  axes  were   essentially  zero.      Stability  of  the   roll  axis   system 
was  not  a   problem  since   proportional  control  was   used   in  the    torque 
generator  and   the   loop   sensitivity  was   relatively  low. 

The   root  Iceus   of  the    linearized   pitch  control   system  was   discussed 
in  Chapter  III   and    shown   in   Fig.    3-3.      The   close   loop  gair,   gave   a 
damping   ratio  of  0.44!:   for  the   two  oscillatory  modes   cf  the    system. 
The   control   system  in   pitch  and  yaw  is  a   fifth   order  system   in  angular 
velocity  as   shown  in  Chapter  III.      The    system  is  also  non   linear,    but 
some   feeling  for  system  characteristics  may   be    obtained   by  comparing 
it  with  a   linear  second   ord >r   system.      Chanter  19  of  Ref.    6  contains 
curves   of  damping  r:  tie  versus   transient   peak  ratio   for  linear  second 
ord^r  systems.      Entering  there   curves  with  the   transient   peak   ratio 
obtained   frc:     the    system  response    giv^s  a   damping  ratio.      A   damped 
perio  ;   of  oscillation  may  also  be   obtain-;   from  the   transient   response. 
An  unaaiu^ad  natural  period  may  be   calculate:!   using  the   damping   ratio. 

The   transient   responses    shown   in  Figs.    5-1  and    5-2  are   not  those 
of  a   second    ,rd^r  system.      Therefore   values   o^"  damring  ratio  and 
undamped  natural   frequency  brsed   on    second   order   res   onse   characteristics 
are  meaningful  only  as  a   rough  measure   of  the    system  characteristics. 
From  Fig.    5-1,    the   yaw  axis   response    indicatas  a  damping  ratio   of  0.48 
and  an  undamped  natural   frequency  of  1.19   cycles   per  second.      From  the 
pitch  axis  response,   Fig.    -"^-2,   a  damping  ratio  cf  0.46  and  an  undamped 
natural    frequency  of  0.95  cycles   per  second  were   obtained.      In  Fig.    5-3 
the    slight   recurvature   of  the    gimbal  angle   curve   indicates  that   the 
roll  control  system  is  rather  heavily  damped.      The   transient   response 


(53) 


indicates   that  the   system  is  dynamically  stable   on  all  axes, 

The   differences  between   inertial,    body,    and   gimbal  angles  are 
plotted   in  Figs.    5-1,    5-2,   and   5-3.      No   interfering  angular  accelerations 
were  applied   to  the   system  during  the   transient   response.      The   existence 
of  these   differences  during  the   transient   response    indicates  that  any 
angular  motion   of  the  missile   ma;/  cause   errors   in  orientation.      The 
absence   of  any  means   in   the    system  of  sensing  these   differences   is 
shown  by  their  constant  value   once   the    gyro  gimbal   angles   have    reached 
their  steady   state   null    position. 

Information  on  the   source   of  the   differences   can   be   obtained   by 
considering  the   orientation  of  the    gyro  axes   relative   to   the  missile. 
To  aid  in  comparison,   the   body  axes  orientation   of  the    gyros  are 
tabulated   below. 


1  Gyro  Axes 

"1 

Rotation  Sensed 

IA 

.y-j\ 

OA 

Roll 

X 

Y 

Z 

' 

Fitch 

Y 

X 

-Z 

Yaw 

Z 

Y 

-X 

J 

Table   5-1 
INTSGRATII.G  GYRO  0RI3KTATIGJJ 

In  the  equations    (2-10)    for  the    integrating  gyro,    the   terms 
introducing  differences   between   body  and   gimbai  angle   are   angular 
accelerations  of  the  missile  about   the   gyro  output  axis  and  the   product 
of  gyro  gimbal  angle   and   angular  velocity  of  the  missile   about  the 
gyro   spin  reference   axis. 

Comparing  the  body  to  gimbal  angle  difference  curves  for  the  roll 
and  yaw  axes  shows  a  marked  similarity  in  the  oscillatory  character  of 
the  differences.      Both  gyros  are   oriented  with  their  spin  reference 

(54) 


axes  along  the   pitch   (Y)   axis  of  the   missile.      This  would    indicate 
that  the  dominant   cause   of  this  difference  was   the    result  of  angular 
velocity  about   the    spin  reference   axis.      The   greater  magnitude   of  the 
error  in  the    roll  axis   plot   can   be   accounted    for  by  the    fact   that,  the 
ginbal  anrle   remains   larger   for  the    roll  axis  than   for  the   yaw  axis. 

The   relatively   small  effect   or;   the   body  to    eimbal  ancle   difference 
resulting  from  angular  accelerations  about  the   output   axis   is   indicated 
in  the   pitch  axis   plot.      The   pitch  axis   gyro  has   its    spin   reference 
axis  along  the    roll  axis.      The    slow   rate    of   change    of   gmbal    angle    for 
the    roll  axis    r~yro   shews   that   the    roll    rate   of   the  rr.issiie    is   low. 
This   indicates   that  the   pitcr.   gyro   is   little   affected   by  angular 
rotation  about   its   spin   reference   axis.      The   output  axis   of  the    pitch 
gyro  is  along  the   yaw   (Z)    axis.      The   curvature    cf  the   Z-axis    gyro 
giir.bal  an^le   curve    indicates   the   magnitude   of   the   angular  acceleration 
of  the  missile   about   that  axis.      It  will   be   ::oted    that  the   peaks   of 
the   oscillations    in   trie   pitch  axis  body  to  gimbal  angle   difference 
curve   occurs  when  the   Z-axis    gimbal   angle    curve    changes  directio:.. 
In  this   transient   response    the   angular  acceleration   of  the    body  about 
a   gyro  output  axis   appears   to   result    in  little   error.      However,    no 
general   statement   seams  warranted    for  the    case    of  arbitrary  missile 
motion. 

The   inertial  to   tody  an;le   difference   curves   are    almost   identical 
for  the    pitch   and  yaw  axes,   while   this  difference    is  almost  nonexistent 
in   the   roll   axis   plot.      This  difference   is   the    result   of  finite    rota- 
tions as  was   shown   in   Chapter  II.      For  the    roll  axis   the   difference 
equation  is: 

(5-D 
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For  th<=  other  axe3  the  difference  is  likewise  the  difference  between 
the  products  of  I ody  angle  an:  angular  velocity  components  associated 
with  the  oth^r  two  axes. 

Comparison  of   Figs,  n-1  and  f>-2  shows  that  the  gimba]  angl^  curves 
are  almost  identical.   This  means  that  at  any  time  during  the  response 
the  body  angle."  and  rotation  rates  are  nearly  identical  for  the  Y  and 
Z-axes.   The  two  terms  Lr.  the  difference  3quation  above  therefore  surr. 
to  almost  zero.   Thi  s  also  makes  the  pitch  and  yaw  axes  differences 
almost  iaentical  with  time  35  shown  in  the  figures. 

The  system  res;  onse  to  an  initial  angular  displacement  could  ha\e 
a  parallel  in  an  actual  missile.   A  similar  event  would  occur  in  the 
case  of  a  second  or  later  stage  of  a  multiple  stage  missile.   In  the 
interval  between  separation  and  ignition  of  the  naxt  stage,  a  transient 
motion  of  the  missile  may  be  expected.   This  transient,  motion  imposes 
the  initial  condition  on  the  control  systei.;  when  the  thru.-t  engine 
fires  and  control  is  commenced.   If  the  angular  velocities  are  small 
and  car.  be  neglected,  this  would  establish  initial  conditions  similar 
to  those  used  for  the  transient  res;  onse.   The  body  v.ould  have  a 
deviation  in  orier.taticn  with  zero  rotation  rates.   For  the  response 
calculated  the  difference  between  inertial  angle  =nd  body  angle  was 
approximately  one  per  cent  of  the  initial  displacement,  for  all  three 
axes.   This  difference  is  the  total  orientation  error  developed  during 
the  response. 

The  information  obtained  from  this  calculation  of  the  transient 
response  of  the  system  is  as  follows: 

1.  The  control  system  is  dynamically  stable. 

2.  Rotations  of  the  missile  rroduce  errors  in  orientation 
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due  to  the  effects  of  three  dimensional  finite  rct-tims  ~n  1  the 
characteristics  of  the  integrating  gyro. 

3.  These  errors  in  orientation  are  net  sensed  by  the 
system  and  exist  in  steady  state. 

4.  For  the  response  computed  the  total  steady  state  error 
in  orientation  was  approximately  one  ;  er  cer.t  of  the  initial  displace- 
ment. 

Periodic  Forcing 

The  systen  response  to  a  periodic  forcing  function  in  interfering 
angular  acceleration  was  obtained.   For  reasons  already  discussed  in 
Chapt-sr  IV  all  of  the  forcing  functions  used  in  this  analysis  were 
composed  of  combinations  of  square  pulses.   For  the  periodic  forcing 
function  the  sequence  of  pulses  shown  in  Fi.-s.  5-4,  5-6,  and  5-8  was 
arbitrarily  selected.   On  a  given  axis  of  the  system  these  pulse: 
are  alternately  positive  and  negative  with  a  period  of  six  hundr-ed 
milliseconds,   iach  pulse  ::as  a  duration  of  one  hundred  mi  1  "J i seconds 
and  a  magnitude  of  four  hundred  milliradians  per  second  squared.   The 
phasing  of  the  pulses  on  the  various  axes  was  arbitrarily  chosen  so 
that  the  yaw  axis  was  twe  thirds  of  a  period  behind  the  pitch  axis 
and  the  roll  a>is  was  two  thirds  of  a  period  behind  the  yaw  axis. 

The  periodic  forcing  function  just  described  is  equivalent  to  a 
linear  vibrational  acceleration  of  the  missile  of  four  feet  per  second 
squared  at  a  distance  of  ten  feet  from  the  center  cf  gravity.   In  the 
absence  of  system  restoring  torque  this  acceleration  acting  for  the 
duration  of  the  pulse  would  produce  a  displacement  of  about  a  quarter 
of  an  inch  at  ten  feet. 

The  response  of  the  system  to  the  periodic  forcing  function  just 
described  is  plotted  in  Figs.  5-4  through  ;;-9  for  all  axes.  These 
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pxots  show  the  forcing  f  unci  ioi  ,  the  gyro  gimbal  angle  ;  nd  the  threw 
angular  differences  for  each  axis  for  the  first  ten  seconds  of  system 
operation.   At  the  end  of  ten  seconds  the  system  had  reached  steady 
state  on  all  axes  so  that  further  computation  was  unnecessary. 

The  plot  of  the  gyre  gimbal  angle  is  essentially  the  same  for 
all  axes.   After  an  initial  transient  pericd  the  gimbal  an^'le  settled 
down  to  a  steady  oscillation  about  its  r  ull  positicn  which  appears  to 
be  almost  sinusoidal.   However,  con  pari sor:  of  these  plots  with  those 
of  Figs.  5-1  through  5-3  shows  that  the  gimbal  angle  response  is 
actually  made  ui  cf  a  continuous  serie:..  of  exponential  curves  similar 
to  those  obtained  for  the  transient  respor.r-e. 

The  behavior  of  the  gimcal  angle  during  the  initial  trar^i-r  t 
period  is  consistent  with  the  results  of  the  transient  analysis 
discussed  previously.   3ach  axis  develops  some  dynamic  error  during 
the  transient  period  which  is  eventually  damped  out.   In  the  case  of 
the  pitch  and  yaw  axes  this  damping  require-  approximately  two  seconds 
while  for  the  roll  axis  approximately  six  seconds  are  required. 

The  effects  of  coupling  among  axes  in  the  system  can  te  seen  in 
the  plots  of  gimbal  angle  for  the  pitcn  and  roll  axes.   Sach  of  these 
axes  is  initially  unforced  and  yet  in  each  case  some  small  gimbal 
angle  develops  immediately  due  to  the  motion  of  the  missile  about  the 
yaw  axis. 

The  angular  difference:  of  rig.  5-4  through  5-9  all  develop  in 
a  fairly  consistent  pattern  in  spite  of  the  vyriation  in  dynamics 
between  axes  of  the  system.   Some  general  discussion  of  the  basic 
patterns  is  warranted  before  considering  the  results  axis  by  axis. 

The  difference  between  body  angle  and  gyro  gimbal  angle  is,  in 
each  case,  an  oscillatory  curve  which  slowly  drifts  away  from  zero. 
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In  steady  state  this  is  a  combii  atioi.  of  a  steady  oscillation  about 
zero  and  a  constant  drift  rate.  Tnis  can  be  explained  in  terms  of 
the  two  error  producing  effects,  discussed  in  connection  with  the 
integrating  gyro. 

In  steady  state  the  component 3  cf  angular  accele ration  of  the 
missile  are  symmetric  about  zero  in  the  present  case.   Therefore  for 
a  full  cycle  the  net  ar.g!-Irir  acceleratioi  cf  the-  missile  about  a 
given  a>is  is  zero.   The  difference  between  body  and  gimbal  angle  may 
oscillate  due  to  the  effect  of  angular  acceleration  about  a  .*yro 
output  axis  but  for  a  full  cycle  the  net  difference  prcduced  by  its 
effect  is  zero.   Any  drift  ir  the  plot,  of  body  to  gimbal  angle 
difference  must  them  te  due  to  the  effect,  of  angular  velocity  of  the 
missile  about  the  gyro  spin  reference  axis. 

The  product  of  gyro  gimbaQ  angle  and  missile  angular  velocity 
about  the  gyro  spin  reference  axis  will  not  in  general  integrate  to 
zero  over  a  full  cycle.   Therefore  over  a  full  cycle  a  net  difference 
between  body  and  gyro  gimbal  angles  wilJ  usually  exist.   In  steady 
state  the  difference  due  to  this  effect  vil3  be  the  same  for  each 
cycle.   This  prcdu  :es  the  linear  inereo.se  upon  which  the  oscillatory 
effects  are  su;  •*  rLr.}  c  .  >  i . 

The  carves  for  the  difference  between  inertial  and  body  angles 
are  characterized  by  a  small  oscillation,  the  mean  of  which  is  linear 
with  a  finite  slope  in  steady  state.   This  linear  divergence  is  to  te 
expected  for  the  case  of  periodic  forcing.   A  net  difference  developed 
during  a  given  cycl^  of  the  forcing  function  would  also  be  developed 
in  any  other  cycle  in  steady  state.   Therefore  there  would  be  a  linear 
accumulation  of  these  differences  with  time  such  as  is  shown  in  the 
curves. 
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The   sum  of  the  differences   just  discussed    is   the   difference 
between  th3   inert  Lai  ar.^1?    and   the    ^yro   gimbal  angle.      This  difference 
is  also  plotted   for  each  axis.      It   represents   the   total  error  between 
the    reference   orientation  and   the   orientation   to  which  the    system  is 
currently  attempting  to  drive   the   missile    iue    to   all  of   the  effects 
under  investigation. 

The  ang  ilar  differences   for  each   axis  at   the   end   of  ten   seconds 
are   given   in   Table    r^-2.      Since   all  .differences    steadily    increase    in 
magnitude   these   are   the  maximum  errors   in  orie-  t;  tion  which  the   system 
experiences   for   the   ten   second    run.      For  those     lifi    rerces  which   are 
oscill'tory  the   numbers   given   are    the   lean  error   at   ter.    swords. 

Table   5-2  also   giv-^s  the    r^tio   of   the   angular   differences   t.     the 
maximum  amplitude    of  t:.-;  missile   angular  oscillation  in   steady   state. 
3ince   the   gyro   sensitivity  from  body   an.-ie    to   gyre   gimbal    sngle    is 
unity,    the   ::.-  -.':-      amplitude   of  gyro   gimbal  angle   oscillation   can  be 
used   in  forming   ratios.      Ko*.;ever,    it   is  rrore  meaningful   to   relate 
errors  in  orientation  to  body   r.r.  :ular  motion.      It   shouli   be   :.    ted   that 
the  magnitude   of  the   body  angle    is  not   that  of  the   gimbal  angle   due 
to   the   errors.      However,    the   amplitude    of  the   body  angle   oscillations 
are   the   same   as   the   amplitude    of  the    gimoal  angle    oscillations. 

The   ten   seconds  of  running  time    for  which   the    results   of   periodic 
forcing  have    been    determined    is    short  compared   to  the   firing  tir.ie   of 
most  missiles.      A   good  apprcximaticn  to  missile   firing  time  would   be 
on;;    hundred   seconds  which  would    require   a   computer   rur.  ten  times  as 
long  as  the   one   that  was  mad?.      However,    once    the   system  has   reached 
steady  state   operation  and    remained   there   long  enough  to  establish 
the   trend  of  the  errors   further  computation  is  unnecessary.      The   error 
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curves,   being  linear,   may  be  extra;   il.'.ted   to  any  desired   value.      In 
the   present  case   this  extrapolation  has  be^-n  carried   out  an  1   the 
resulting  errors   :n   orient  tion  are    shov/n  in  Table    F— 2  alor.g  with 
the   values   for  ten   seconds.      The   values   of  the    differences  at   one 
hundred   second:-    ire    given   in   both   radians  and  degrees   for  convenience. 

The   following  infonnati    r.   was   obtained   f r<  rr.   the    response   of  the 
system,   to  a   periodic   interfering  angular  a  :celeraticn: 

1.  Errors   ir.   orient-tier,  developed    due    to   the   effects   of 
thre3  dimensional   finit?   rotations   and   the   characteristics   of  the 
integrating  gyro. 

2.  These  errors  v.jr''  sys4.  ?.T.atic ,  an'J  in  steady  state,  their 
mean   in:re^-sed   linearl;    with   time. 

3.  The   total   error   in  orientation  develop;.1   at   the   end   of 
10  seconds  v;as  0.'A1   millircdians   in   roll,   1.05  millir;  Jians   in   pitch, 
and   1.4?  milliradians   in  ya-.  .      These   art-    respectively   34b,    27^,    and 

40'o  of  the   arr  litude   cf  missile    angular  ;.;oticn  about   the   axis   in  question. 

4.  Jince  the  increase  in  error  was  linear,  extrapolation  tc 
100  seconds  was  _ossible.  This  increases  the  values  above  by  a  factor 
of  ten. 

Random  Forcing 

In  order  to  simulate    interfering  angular  accelerations  as 
realistically  as   possible   witbra   the   limitations   already  imposed  by 
computing  time,   the   response   of  the   system  tc  a    series   of  random 
pulses  was  determined.      Since   no   steady  state   could   be   expected   for 
such  a    fcrcing  function  the   problem  was   run  for  forty  seconds  of 
real  time.      This   required  approximately  ten  hours   of  computing  time. 

The   random  forcing  function  was  made   random  with   resrect  to 
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three  of  th?  pulse  cha ract-;rLstic3.   The  sign  >:"   Ih-  pulse,  Its 
magnitude,  and  its  duration  w;h  independently  chosen  in  a  random 
fashion.   This  was  dona  automatically  in  the  computer.  The  sign  of 
the  pulse  '..as  chosen  first.   The  mechaidsm  :  y  which  this  was  done  was 
such  that  positive  and  negative  signs  each  had  a  probability  of 
occurrence  of  one  half  in  a  riven  selection.   The  magnitude  of  the 
pulses  was  limited  to  the  values  from  0  to  90  milliradians  ner  second 
squared,  in  steps  of  ten  milliradians  per  second  squared.   ;Cach 
value  had  a  probability  of  occurrence  of  cne  tenth  in  a  given  select- 
ion.  The  duration  of  the  pulses  was  selects  at  random  in  the  range 
from  5  to  500  milliseconds  witn  selection  being  llmitel  to  multiples 
of  five  milliseconds.  The  forcing  function  obtained  in  this  manner 
is  shown  for  each  axis  in  rigs.  5-10  through  5-12.   It  should  be  noted 
that  there  is  n^  correlation  between  the  forcing  functions  for  the 
various  axes  since  selection  or.  each  axis  was  completly  independent 
of  the  other  two. 

The  response  of  the  system  to  the  random  forcing  function  just 
described  is  shown  in  Figs.  5-10  through  r-13. 

The  system  response  as  shown  by  gyre  gimbal  angle  is  about  what 
could  be  expected  in  view  of  the  previous  results.   The  gimbal  angle 
oscillates  randomly  in  essentially  the  same  manner  as  the  forcing 
function.   There  is  no  steady  stat?  condition  apparent  in  the  response. 
None  was  expected  since  the  forcing  function  is  random.   In  general 
the  behavior  of  the  gimbal  angle  is  consistent  with  the  rrevious 
observations  on  system  dynamic  stability. 

The  angular  differences  have  the  same  general  character  as  those 
observed  in  the  periodic  case.   The  body  to  gimbal  angle  difference 
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curve  is  oscillatory,  while  the  inertia!  to  body  angle  difference 
is  relatively  snooth.   However,  there  is  no  systematic  divergence 
of  these  angular  differences  for  this  run.   In  this  respect  the 
response  was  different  from  the  results  for  periodic  forcing. 

3ir.ce  the  response  to  random  forcing  is  essentially  statistical 
in  nature  generalized  statements  Ccn  not  be  made  on  the  b3sis  of  a 
single  run.   The  angular  difference  curves  shown  are  just  one  set  of 
an  infinite  number  that  could  occur  as  a  result  of  this  type  of 
forcing.   The  fact  that  the  angular  differences  obtained  during  one 
run  do  not  diverge,  does  not  imply  that  they  would  not  diverge  in 
another  run.   however,  the  single  run  shown  does  demonstrate  that 
errors  in  orientation  do  exist  ever,  with  random  forcing.   The 
interfering  angular  acceleration  inputs  nee!  •  ot  be  periodic  for 
these  orientational  errors  to  develop;  And  these  are  the  or.ly  conc- 
lusions th  t  C3n  reasonably  be  drawn  fro::,  the  res-ilts  of  a  single 
run.  More  runs  would  permit  a  statistical  approach  to  the  problem. 
More  generalized  results  could  then  be  obtained. 

A  run  was  made  to  show  that  orientational  errors  will  exist 
regardless  of  control  system  dynamics.   The  system  dynamics  were 
altered  by  assuming  proportional  response  in  the  integrating  gyros 
and  the  torque  generating  systems.   The  reduced  equations  are  shov.-n 
in  Appendix  B. 

This  simplified  system  was  forced  by  the  same  periodic  forcing 
function  used  previously.  The  response  to  this  forcing  is  shown  in 
Figs.  5-14  through  5-16. 

It  will  be  not^d  that  the  plots  of  this  response  show  missile 
body  angle  as  the  reference  variable.   For  this  simplified  system 
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the  gyro  angle  is  identical  to  the  body  angle.   Therefore  the  only 
difference  that  can  occur  is  that  between  inertial  and  body  angles. 
The  angular  difference  curves  from  inertial  to  body  angle  ir.dicate 
the  error  in  orientation  develops  in  the  same  manner  as  for  the 
complete  system.   The  error  is  slightly  less  in  this  case  3ince 
the  control  system  response  is  faster.   Therefore  the  missile  will 
experience  smaller  angular  deviations  and  lower  angular  velocities 
for  the  same  forcing  function.   This  accounts  for  the  reduction  in 
the  net  error  in  orientation  due  to  finite  rotation  effects.   In 
any  case  the  error  does  exist  regardless  cf  control  system  dynamics. 

The  results  obtained  in  this  investigation  prove  the  existence 
of  errors  in  orientation  due  to  three  dimensional  finite  rotations 
and  integrating  gyro  characteristics.   Some  information  vss  obtained 
for  the  cases  of  periodic  forcing,  random  forcing,  and  transient 
response  from  an  initial  angular  displacement. 

For  the  case  of  periodic  forcing,  further  investigation  is 
needed  to  determine  the  effects  of  the  frequency  and  magnitude  of  the 
forcing  function.   In  the  case  of  random  forcing  more  runs  will  be 
required  to  permit  a  statical  evaluation  of  the  errors  developed. 
If  very  many  runs  are  to  be  made,  a  faster  mehtod  of  computation  is 
highly  desirable.   This  can  be  accomplished  through  the  use  of  a 
faster  computer,  or  perhaps  by  changing  the  method  of  computation. 
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CHAPTER  VI 

ca;ciU3ic:;s  a:.:  rsccm-skdations 

The  preceding  analysis  has  shown  that  errors  ir  orientation 
may  develop  in  the  geometrical  stabilization  >ysteir  using  body  fixed 
^ros.  These  errors  are  iistinct  from  the  levial  '  n  In  orientation 
which  such  a  systerr  would  ex.-erience  due  to  system  dynamic  lags.  It 
was  demonstrated  that  under  certain  conditions  the  errors  in  orient- 
ation systematically  increase  v.-ith  time.  AIs;  ,  it  was  shown  that 
errors  in  orientation  develop  regardless  of  system  dynamics, 

Zrrors  in  orient- ticn  developed  in  the  absence  cf  interfering 
angular  vibrations.   In  restoring  from  an  initial  anrjlar  displacement 
the  system  developed  an  error  of  about  one  per  cent  of  the  initial 
displacement.   This  error  remained  when  the  system  had  come  to  rest. 
Those  errors  developed  due  tc  the  coupling  of  missile  angular  motion 
between  axes  of  the  system. 

'..hen  the  system  was  forced  with  a  periodic  interfering  angular 
acceleration  or.  all  xs,  there  was  a  systematic  increase  in  orientation 
errors  with  time.   After  ten  seconds  the  total  error  in  orientation  in 
yaw  was  forty  per  cent  of  the  maximum  amplitude  of  the  angular  oscillation 
of  the  missile  in  steady  state.   Corresponding  figures  fcr  the  pitch 
an  J  roll  ax.es  were  twenty  seven  .  er  cent  and  thirty  four  per  cent 
respectively.   These  errors  increase  linearly  with  time.   Therefore  at 
one  hundred  seconds  the  error  in  yaw  becomes  four  hundred  per  cent  of 


the   maximum  amplitude   of  the   angular  oscillation  of  the   rrdssile. 

The   same   kind   of  errors  in  orient.- tion  developed  when  the 
interfering  angular  sc  cell  rat  ions  were   random.      However,    due   to   the 
statistical  nature   of  the    forcing  function  n:-   conclusions   could   be 
made   on  the    basis  of  a   single    run   concerning  the   manner  in  which  the 
errors  will   accumul^t^    with   time. 

rurth*r   Investigation  cf  th^   :  roblem   is   still   required.      The 
investigation  of  periodic   forcing  should   be   exten'ed    tc    determine 
the   effects  of  the    frequence  ar.d   magnitude   cf  the    interfering  anguia: 
arc^ler  tions.      Jr..  ugh    results    should   be   obtained   using  a    rand    rc 
forcing  function  to   penr.it  a    statistical    analysis   of  the  errors. 
These   additional    investigations  v.'ouli    be    greatly   facilitated    if  a 
significant   reduction   in    co-out  in  ■  tin  •    :ould   be    achieved. 
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API £KDIX  A 
TH3  TRANSFER  OF  VICTOR  CCl-.FCr-E.'TS  EETWESN  ROTATED  AX33  SITS 

Fig.  A-l  shows  the  inertial  and  body  axes  sets  de scribed  in  Chapter 
II  and  ar.  arbitrary  vector  V  which  is  fixed  in  the  body  axes  set.   The 
:  roblem  is  to  jxpress  the  components  of  the  •/•  tor  in  the  inertial 

sat  when  the  \   xh  axes  components  are  known  and  the  body  axe? 
have  assumed  an  arbitrary  orientation.   In  general,  an  arbitrary 
orientation  of  the  body  axes  relative  t<  the  inertial  a  es  can  be 
duplicated  by,   t  most,  three  successive  finite  rotations.   In  the 

lysis  be  lew  three,  successive  rotations  c-f  the  body,  one   about  each 
of  the  inertial  axes,  are  assumed.   After  each  rotation  the  components 
cf  the  vector  V  in  inertial  a.es  are  derived.   After  three  rotations 
the  results  for  the  individual  rotations  are  combined  to  give  the 
components  of  the  vector  in  inertial  axes  in  terms  cf  the  three 
angles  of  rotation  and  the  components  in  body  axes.   The  resulting 
transformation  equation  is  not  unique.   There  is  a  different  trans- 
formation for  each  possible  order  of  rotations.   However,  a  simplification 
is  made  at  the  end  which  gets  around  this  difficulty.   In  the  derivation 
below  the  order  of  rotations  is  about  the  •'• ,  T ,  and  Z-axes  in  th.^t 
ordar. 

Let  the  body  axes  of  Fig.  A-l  be  rotated  about  the  X  inertial 
axis  through  the  angle  A„  as  shown  in  Fig.  A-2.   From  the  figure  it 

(80) 


can  bo  seen  that  the  components  of  th3  vector  V  in  inertial  axes 

can  now  bo  written: 

(A-l) 

V  =  V 
*I    *B 

V  =  V  cos  A    -  V   sir.  A 
yI    yB     *I     ZB     *I 

V  «  V   sir.  A    +  V   cos  A 
2I    *B     XI     Z3     XI 

There  is  considerable  simplification  -iffcr-:\>;i  by  writing  the  set  of 

equations  (A-l)  ir.  matrix  form. 


(A-2) 


v   "!     'l 


*I 


V  =     0 

yI 


0  0 


V 

xi       *i        yB  : 


!  0  5A 


x^ 


V 


■: 


B 


Where   the   notation  5A      t   and   CA      refeis      to  the    sine   and   cosine   of 
the  angle   of  rotati  >n  about  the   X-axis. 

Nov;  let  a   new   set   of  body  axes    (8   }    he   defined    such   that   they 
are   coincident  with   the    inertia!   axe:;   of  Fig.    A-2.      I.'ote    that   the 
transformation   (A-2)   always   applies   between   the    sets   of    -oiy  axes 
(B)    and    (B').      That   is: 
(A-3) 


V 


0  0  v      ! 

*B  i 


=     0  CA        -3A  V 

*I  *I  yBi 


ZBi 


0  3A  CA  V 

X-  X_  2 


B 


J     L_ 

If  the   body  is  now  rotated  about  the   Y  inertia  1  axis   through 
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A   then  from  Fig.  A-3  i>.  can  be  seen  that: 


(A-4) 


1 

L_    — 


CA 


3A 


7-1 


-3A 


CA 


V 


and  definir.g  the   body  axes   set   (5,!)    in  exactly  the    same  manner  as 

(3«)    gives: 

(A-5) 


— 
I    V  ' 


V 
*B 

V  i 
^8 


V" 
Z3 


,    CA 


:-i 


SA 


*B 

B 


V 


-OA 


0  CA 


'B 


..0 


(A-6) 


tating  nov:  about  Z_  through  A        gives: 

1  Zj 


r 


*I 


CA 
z-, 


-5A 


CA 


XB 


V" 


0 


B 


V" 
ZB    ! 


Replacing  the  matrices   above   by  symbols   of  the    form   (V-r)   allows 


equation 

(A-6) 

to 

be   v;r: 

Ltter 

(A-6) 

(V 

=: 

<Az> 

(V  ' 
v   B 

) 

where   the  meaning  of  the   notation   becomes  evident   by   comparing  the 
^wo  forms  of  the   equation   (A-6). 

If  (A-3)    is   substituted   into   (A-.1:)   and   this  equation  in  turn   is 
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substituted  into  (A-£)  the  result  is: 


(A-7) 


(ty  -  (A2)  <Ay)  (Ax)  (VB) 


If  the  matrix  multiplication  (A  )  (A  )  (A  )  is  carried  out  the  resulting 
equation  for  the  inertia]   components  in  terms  of  body  cor.;^onp:.t s 
becomes: 


(A-*) 
l~V 


*I 


CA  CA 
2I  yI 

JA  CA 
ZI  yI 


*,A     3A  3A   -SA  1h 


'l     yI  *I 


ZI  >:T 


CA   SA   CA   +3A   3A   '  'v 
ZI  yI   XI    ZI   *I 


V 


L 


3A   3A   +CA   CA 

:I    ZI   ^1 


yI   *I  'Zt    Zt 


*     r»  * 

:I  yl  *i 


ZI  xi 


yB 


CA 


CA   CA 


-I 


V 

L  !i 


xi 

If  the  rotations  about  the  inertia]  axes  ar  kept  s-all,  so  that 
the  product  of  the  sin  3  of  two  of  the  angles  can  be  neglectei  this  can 
be  written: 


(A-9) 


*I 


[~i 


fv  1 


1 


-A 


XI 


-A 


*I 


V 

ZR 


It  will  be  noted  that  the  restriction  imposed  on  the  size  of  the 
angles  in  (A-°)  is  more  severe  than  that  imposed  by  the  usual  small 
angle  assumption.   Therefore  replacing  the  sine  of  the  angle  by  the 
angle  and  t,he  cosine  by  unity  is  justified. 

It  wa3   noted  above  that  the  transformation  given  by  (A-")  is  not 
unique.   In  fact  there  are  six  different  matrices  similar  to  the  one 
in  (A-3)  corresponding  to  the  six  orders  of  rotation  possible  about 
the  inertial  axes.   However,  under  the  small  angle  assumption  made 


(83) 


these  all   reduce   to   (A-?).      The   equation   (A-9)    and    (2-2)    are    identical. 


(S4) 


^z 


8 


VH 


*Ii*S 


Initial    0*iB*Tf\TiON    of     Boor  A*cs,X«£«Tial  Ax£Sj  a*o  ^cro/f   V 
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**.* 


O*  IE  NT AT/o*)     ArrCR       Rqtat/oa?   A  Boot     X    IvERTiAL    Ax'5 

F,g.   A-ji 


(36) 


(s 


Orientation  Aft  BR  Rotation   Asout     Y Jnbrtial   Axis 


(87) 


f//v*L   Oricntatioh  After   Rotation  As out  Z.  £*£/?Til  Axes 

F.%  A-f 
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AFFENDIX  B 

SYSTEM  PERFORMANCE  STATIONS 

The   performance  equations   for  the    system   cf  Fig.    2-1  are   developed 
belovv   in  the   order  of  their  discus:. ion   in   Chapter  III. 

The  missile    is  taken  as   a   rigid    body.      Since   translational  motion 
is  of  no  interest   in  this   analysis,    the   center  cf  mass   of  the  missile 
is  assumed    fixed    in   inertial    space.      The      is~ile  equations   of  motion 
are    then  obtained   by   summing  moments  about   the  missile   certe r  of 
gravity.      Expressing  components   in  body  axes   gives  Zuler's  equations: 
(B-l) 

W     +   (I     -I   )  W  W     =         K 

x  z       y       z  y  x 

W     +   (I     -I    )  W  W     =         M 
y  x       z        x  z  y 

W     +   (I     -I    )   W  V     =  M 

z  y       x       y  x  z 

where  in  the  present  case  the  sum  of  the  applied  moments  on  the  right 

above  includes  the  restoring  torque  applied  to  the  missile  by  the 

control  system  and  the  interfering  torques  due  to  missile  angular 

vibration.   Thus: 

(B-2) 

M  =  K   +  !•:,.  .x 
x    m     Cint) 
x        x 

y    m     (int; 

z    m     Cint; 
z        z 
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The  ren'.aind'T  of  the  system  performance  equ.otion:3  are  concerned 

with  the  determination  of  the  restoring  toraue  (K  )  as  a  function  of 

*     ra 

missile  motion. 

In  the  missile  configuration  being  considered  restoring  torques 
for  the  pitch  and  yaw  axes  (y  and  z  axes)  are  produced  by  rotating 
the  thrust  line  of  the  missile  engine  away  from  the  nissile  center 
of  gravity.   Restoring  torque  is  roll  is  produced  by  two  nozzles  on 
the  circumference  of  the  missile  at  the  opposite  ends  of  a  diameter. 
The  components  of  restoring  torque  about  the  pitch  and  yaw  axes  may 
be  written: 
(B-3) 

M    =  T  1  A 

m         m 

y      y 

K   =  T  1  A 

m         m 
z         z 

Where:     T  ■»  the  thrust  of  the  missile  motor,  assumed  constant. 

1  ■  distance  from  the  missile  center  of  gravity  to  the 

motor  pivot  point. 

A   =  the  angle  through  which  the  rotor  is  rotated  about 
m         o        o 

the  pitch  or  yaw  axis  (Assumed  small) 

Similarly  the  restoring  torque  in  roll  becomes: 

(B-4) 

K   •=  T,   v  d 
m^    (  re  ) 

'/."here:     1*      x  =  the  thrust  of  the  roll  engine, 
d     =  the  diameter  of  the  missile. 
The  system  which  positions  the  missile  motor  in  pitch  and  yaw 
is  assumed  to  be  a  second  order  system  which  receives  gyro  signals 
as  input  and  positions  the  missile  motor  about  the  appropriate  body 
axis. 
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The   performance  equations  may  be   written: 
(B-5) 


P         2 


2(DR) 


3 


V."  — .,     " —     p  +   ll  A  ■  S        /  k\^f.      \ 

Pes  npcs  J      y  igu'  y 

pcs(e         :A)    (rru) 
r         rgu  y 

2 
n  — ,,     •>  -~  +   II  A  ,\e/. 


3 


yes  V.  I    m  ycs(e,       ;A)    (igu) 

yes 

ycs(e        ;A)    (rgu) 


..;.ere: 


(Da)  =  7  itch   cc.-.trol   system   damping   ratio. 

pC  5 

V.'  =  }  itch   control   system  natural   frequency. 

*i;cs 

3   ,  ,  >  =  '•  itch  control  svstem  sensitivity  from 

ccs(e.   ;A) 

igu 

integrating  gyro  signal  to  motor  angle. 


a 


ocs(e        :A)    =   Fitch  control   s  -sterr-.  sensitivity   from 
rru* 


r-^ta   gyro   signal  to  motor  angle. 
Similar   definitions   hold    for  the  yaw  ccr.trcl    system. 

.loll  engine    thrust   is  assumed   to  be   directly   proportional  to 
gyro   signals.      Tnat   is: 
(3-6) 

(re)  rcs(e.       :*)    igu  rcs(e        ;7;    rgu 

igu*  b  x  rgu  °  >: 

Scuatio:  s   (3-3)    and   B-5  may     e    combined    to  give   the    restoring  torque 

in   pitch  and  yaw  as  a   function  of  gyro  signal.      2ouati*>r.s(?-4)    =nd 

(E— 6)   may  be   similarly  combined.      Then: 

(B-7) 

'in     *     (rcs)(e.      :K)    igu      '      (rcs)(e        ;M)    rgu 
x  igu*  &  x  rgu*    '      °  x 

(91) 


p 
w 

n 

_       pes 


yc3 


2(DR)  . 

J£S    p       ♦    1 

v: 

n 
pes 


2(Dh)  . 

♦  ^cs   p     +   1 

W 

n 

yes 


M 


m 


a         ^ 


pcs(e.      ;!•')    igu 
igu'  b  i 


pcs(e        ;TJ    rgu 
r         rgu*  &  y 


/  \  e 

ycs(e.       ;K)    igu 


ycs(,e        ;>.>    rgu 
J  rgu'  6  z 


Whe  re : 


3(rcs)(e,       ;!•:)    =  d  S(rcs)  (e,       :T) 


igu 


igu 


3(yca)(eiguil.;)   -  T  1  S(yc3)(e;  j 


IgU' 

The   remaining  sensitivities  have    b-?en   simiis ■  rly   red-fine. 1. 

In  writing  equations  to  describe   th*   gyro  units,   the   gyre 
uncertainties  are   neglected   due   tc   the   relatively   short   systerc 
opera  ting  time.      The   rate    gyro  units   &re   assured   tc  be   proportional 
mechanisms   so  that: 
(B-8) 

rgu^  rgu(/.;e)      x 

rgu^  rgu(l.'je)     y 


e  -     S        ,,,     x   W 

rgu  rgu(Y.;eJ      z 

where   the   sensitivj t.v   i3  trie    rate   gyro   sensitivity  from  input  a>:t> 
angular  velocity   to   outrut   signal. 

The  performance  equations   for  the   integrating  gyro  is  developed 
in  Derivation  Summary    2  of  Reference    2  .      In  terms  of  gyro  input, 
output  and   spin  reference   a>es,    the   torque    summation   on   the    gyro 
gimbal   is: 


(92) 


(B-V) 


^(gim)   A(gim)  +  A(I-B)CA  ~  A(vd)(A;K)  A(gin.) 

•  • 

♦H  A/T  DwT1>co3  kf        \  -  H  A/T  nw^,vsin  k,    .  \  -  0 
(I-B)(IA)     (gut)      (I-B)(oRA)     (gim) 

The  orientation  of  the  gyro3  with  respect  to  body  axes  is  given 

in  the  following  Table: 


Gyro 

In 

>ut   Axis 

opm   Reference 

Axis 

Out} 

ut   /xis 

Roll 

X 

Y 

Z 

Pitch 

Y 

X 

-z 

Yaw 

Z                                          Y 

-x 

Table   5-1 

3y  :••:■  Crie 

sntaticns 

where   the  entries   in  t'..e   t-jble    refer  to  the  missile   body  axes. 

Trie   relationship  between  gyro   gur.bal  angle   and    gyro  outrut 
signa]    Ls: 
(B-10) 

igu  '"      (igu)(A        ;e)      gim 

(B-9)   -^nd    (-—10)   xay  be   combined   tc   give   the   equation   for   gyro 

outj.ut   signal,      Equations  ir.ay  be   written   for  each  axi?   by   replacing 

the    subscripts   refering  tc   gyro  axes  by  the  appropriate   body  axis 

subscripts   frc.T    the   Ta'--le    (B-l).      Then,   after   rearranging   (3-9) 

becorr.^s: 

(B-ll) 

••  • 

If  .,    \efj  ___>      +  3/..jN/t..  .,,8,..      x      +HW.e 


(gimrdeu:.  "(vd)0V;K)°(igu). 


y    (igu). 


=   H  S,„ 


v;    -  i 


(Sg)(A;e)     y  "      (gim)°(sg)(k;e)      z 
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I(gln)e(igu)     +  S(vd)(V:.M)e(igu)     *  H  Wx  e(iKu) 

9 

"  H  3(Sg)(A;e)Wy  +   I(gi3i)   S(Sg)(A;e)  Wz 

•  •  • 

(£i.T.)e(igu)  (vd)(V/jK)e(igu)  "y  e(igu) 

• 

=  H  S(Sg)(A;e)Wz  +  I(gim)   S(Sg)(A;e)   Wx 

The  equations    (D-l),    (B-2),    (B-'7)    and    (B-ll)    completely   :••    cribe   the 
jerforu.ance   of  the   proposed   system.      These    can  be    simplified   by 
redefining  certain   of  tr:e    constants. 
Define : 


(B-L2) 


I      -  I 

(IR)     --*- 2 

x  I 

X 

i    -  r 

(IR)      =      ZT        * 
y  I 

y 
i    -  i 

(IR) 


z  I 

z 


/-.  .L3C   xet: 


(int) 


in  •« 
x       ,                         x  , 

I  "      (int)  I        =     ra 

x  x   ;        x  x 

K(int) 

-*-  A,. 


>: 

.i. 

•  • 

V 

_"" 

=  «. 

1 

m 

y 

y 

I  ■  "(int) 


int  n 
_z        .                       _z  . 

I_  :      (int)  I        =      (int) 


z  z    ;        z  z 

The   gyro  angular  momentum  and   the   viscous  damper   sensitivity 
are   equal  for  the   gyro   beir.g  used    so  that: 

(B"14)  I    .  I    . 

(CT).        -^  --gffl 

^       S(vd)(W;M)  H 
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Substitution   of   (b-12)    (B-13)    (B-14)    into   the    system  equations  hag 

the   effect  of  lumping  the    system  sensitivity  for  each  axis   in   the 

torque    genera  ting  systen.  equations,   except   for  the    signal   generator 

sensitivity.      If  the    gyro  equations  are  expressed    in   terms  cf  gyro 

gimbal  angle   instead   of  voltr^e   this   sensitivity  is  also  carried 

up   into  the   torque    generating  system.      The    perform   nee   equations 

for  the   systems  then  be ecus: 

(B-15) 

.  . .                   •  • 

W     -   (IR)   '..  W  =  A, .    .  x      +  A 

x              'x  y  z  ^lnt^x         mx 


W     -  (IR)  v:  V.      =  A,.    4>      + 
y  y  x  z  (int.) 


•  • 

A 

m 


W     -   (IR)   '..'  V.'     -   A, .    ,  x      + 
z  z  x  y  v,  int) 


m 


...  +     n  *'        V 

m     "       (rcs)(A    .      A)      gim  rcs(?.:;A)     'x 


!(DR) 


d  :  s 


P  -  1 


',^CS 


1C3 


rr.  pcs(,n    .     :A)      gun 

y         *         gun'  -     y 


2(DR). 


yes 


'CS 


+   S 


>pca(V:;A)    "y 


yes  (A    .     :A)       .-;i.x 
J       v    gun'    '       °      z 


+   S 


pcs(V.';A)    "z 


(CT).        A    .        +  A    . 

igu     gimx  gunx 


;  A   .        =  v;     -  (CT). 
y   ?un  x  1  ?u     z 


(CT).        A    .        +  A    .         +  V;  A    .        =  W     +   (CT). 

'xgu     gxm  gun  x  guuy         y  'igu     z 


(CT).        A    .        +  A    .        +  W  A    .        =  W     +   (CT). 

y  gin  z  igu     x 


igu     glm2         gln-.z 


gim  I  . 
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~  (rcs)(e        :M)      (rgu)(W;e) 

S(rcS)(v;;A) j& 

x 

and  identical  definitions  for  the    sensitivities  held   for  the    remaining 

axes  of  the    system. 

The  equations   (B-15)   are   the    set  that  wa.     solved   on  the   digital 
computer  in  the   present  analysis.      The   values   chosen   for  the   various 
constants   in  these   equations  are   given    L:    Tatie   3-1. 
Reduced  Dynamics  Case 

The   sir.  piified    system  which  was   Jiscossei    in  Chapter  ITI   can  be 
obtained   from  the   equations   (lr-15)    b\    eliminating  the   dynamics   and 
the   non-linearities   fro::,  the    last   six  equations.      These   equations 
can  then  be    coo.bir.sd   with  the   missile   equations   of  rr.ction  to   give: 
(B-16) 

\  ~   (lR)xV^  =  AUnt)      +   S(rcs)(A„^;A)      A(gi:  ) 

X  &i..  X 

(rcs)('.;;A)    "x 


W     -   (IR)   \i  V     =.  A/,    .v      +  S/         w.  \\      A/    .    x 

y  'y  x  z  Unt)  (pcs)(A    .    jA)        ( gim) 

j  &•»-'•  y 

(pcs)('..';A)    "y 


V.:     -   (IR)   Vj  VJ      =   A/.     .  n       +   3/         w,  *\      A,    .    v 

z  'z  x  y         (int)z  (ycs)(n        ;^)        (gim)z 

(pcs)('..';A)    "z 

V-Tiere:  A,    .    >      ■  V* 

(gxm;^  x 

. 
(gim)y         y 

(gim)      =   "z 
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The    set  of  equatior.5   (B-16)    waj  the   set   solv«:>l   for  the    system  with 
simplified  dynamics.      The   constants  here   identical  with  those   for 
the    cor.plete    system  which  are    given   in  Table   3-1. 
Single   Axis,    Linearized   Stability  Analysis 

The   question  of  systerr    stability  was   a;  proached   by  making  a 
root  locus  plot   of  the   linearized   performance   equations   for   the 
yaw  axis.      The  yaw  axis  equations  may  be   written: 
(B-I7) 


Wz  -    .92?  v;x„..  ♦   A(.nt)        *     At 
J  z  z 


p 
"525 


-~     +   1    I    A 
<cT  I      m 


J, 


ycs(A    .    A) |       (gin) 


_     ycs(.-.;A) 

Sycs(A    .   A) 
gun 


r 


p 

!l33 


+•       D       + 


gun 


=    V,'  4    -==- 

z         333 


.^r^r?  numbers  h   ve    bean   substituted    from  Table   3-1   for  all  except   the 
sensitivities.      These   may  tn   combined   t:>  give   a   single   expression 
provided  th,-.t  the   character  of  the   t.on-lineorities   is   preserved. 
That   is,   the   order  of   factors   in   the   non-linear  expression  must  be 
maintained.     Then  by  substituting  from  the   first  and   thin!   equation 
into  the   second   and   their  results: 
(  B-1S) 


333    +  p  +  V 


*P2 
|o25 


~.     ♦   1 


(ocs)(A    .    ;A)    i 


333 


W     -   .925-:  W     -  Ax.    .x 

z  x  y  (mt). 


1  ♦  .1985^  ♦  p  +  vj      V.-!- 


(9') 


where   the    rati 


(B-19) 


VcsQ.;A) 

Vcs(A    .    ;A) 

'  Rim 


-      .1925 


is  arbitrarily  chos?n 

If  it  is  now  assumed  that  there  is  no  motion  In  roll  ?nd  pitch 
then* 
(B-20) 

w      =  v;      =o 


The  equation   (3-]?)    then   becomes   linear   ~nd   the    performance    function 

for  tne   yaw  axis   rr.ay   be    written: 

(3-21) 


i?W. 


int;"x 


2 
+   P     To*  +  -fe  +   ii 


333 


J 


3§3  ♦  1 


£_  +  £_  + 
625       25 


+   i> 


:s(Ag.m;A)L_333 


.1QS5P  +  1 


and  the    characteristic   equation  is: 
(B-22) 

,       Sycs(.19'*5  x  .003  p2  +   .1?*."  p  +  l)  n 

P^r.oo3p  + 1  pf  +  2_ 

L  -1625       25        "_] 

=   1  -   (IF)  .  =      0 

open   loop 

The  complex  frequency  plot  of  the  open  Icoi  performance  fur.cticn 

(PF)  ,  as  a  function  of  3   /      v  is  given  in  Fig.  3-3.  The 

*  gim' 

plots  and  zeros  of  this  plot  are  given  below. 


(9*) 


Poles  Source 

P-.    «   -333.3  Integrating  Gyro 

Pp,F„   =  0  Integrating  Gyro  and  Kissi.le    equation 

P4,Ft-  =   -12.^     j  21.65  Yaw  Torcue    Generating  System 

Zeros 
Z     =   -32s?. 9  Hate    Feedback 

Z2  «  5.1  -late   Feedback 

As   the    open   loo;-   sensitivity   is    increased    from   zero   it  will 
be   observed   that   the    pole    pair  at   the      rgin    .  .  ve   to  the   left  while 
the   pole   pair  due    to  the    tor:.--    generating    >ystenr    nave    to   the 

right.      The   choice    of   the   open   loo:     sensitivitv  3        , ,  .  >    was 

*   ycs(Ag.m;A; 

made  so  as  to  give  both  of  these  oscillatory  modes  the  same  dam:  ir;g 
r-tio.   This  rrear.s  that  both  cries  rr.ust  lie  on  the  same  radial 
line  fro:;;  the  origin.  By   trial  and  error  this  radial  was  deter  ir.ed 
to  ^a  the  line  for  a  damping  rati:  of  D.44.   When  both  pole  pairs 
lie  on  this  line  the  open  loop  sensitivity  was  deter?  ined  as: 
(3-23) 

Sycs(A3in;;A)   =  — ° 

which  is  the  value  shown  in  7-?bie  3-1. 
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AFFiCNDIX  C 


Cr  y     i-   ■  .  -    -,r      -  •      Mi 

.   .   u  i  ."jit    I  ...     rru  . 


Trie   foil  'A'ing  program  is  the    complete    program  us«d   for   the 

system,   with   random   forcing  functioi  .      The   program    is   coded   us 

(14U15) 

the  FLAP  systen.    ;v    .   Other  programs  used  are  simply  modifica- 
tions of  this  basic  program. 
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